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Abstract 



Using a bidifferential graded algebra approach to "integrable" partial differential or difference 
equations, a unified treatment of continuous, semi-discrete (Ablowitz-Ladik) and fully discrete 
matrix NLS systems is presented. These equations originate from a universal equation within 
this framework, by specifying a representation of the bidifferential graded algebra and imposing 
a reduction. By application of a general result, corresponding families of exact solutions are 
obtained that in particular comprise the matrix soliton solutions in the focusing NLS case. The 
fyr") ' solutions are parametrised in terms of constant matrix data subject to a Sylvester equation 

fT^ , (which previously appeared as a rank condition in the integrable systems literature). These 

data exhibit a certain redundancy, which we diminish to a large extent. 

More precisely, we first consider more general AKNS-type systems from which two different 
matrix NLS systems emerge via reductions. In the continuous case, the familiar Hermitian 
, conjugation reduction leads to a continuous matrix (including vector) NLS equation, but it is 

well-known that this does not work as well in the discrete cases. On the other hand there is a 
complex conjugation reduction, which apparently has not been studied previously. It leads to 
square matrix NLS systems, but works in all three cases (continuous, semi- and fully-discrete). A 
^> , large part of this work is devoted to an exploration of the corresponding solutions, in particular 

■ regularity and asymptotic behaviour of matrix soliton solutions. 
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1 Introduction 



Using the framework of bidifferential calculi (or bidifferential graded algebras), a unification of 
integrability aspects and solution generating techniques for a wide class of integrable models has 
been achieved [1] (see also the references cited therein). In particular, there is a fairly simple and 
universal method to construct large families of exact solutions from solutions of a linear system of 
equations. The formalism treats continuous and discrete systems on an equal level. This suggests 
a corresponding study of the well-known nonlinear Schrodinger (NLS) equation (see e.g. [2,3]), its 
semi-discretisation, the integrable discrete NLS or Ablowitz-Ladik (AL) equation [4-7], and a full 
discretisation [5] . The literature on mathematical aspects of the continuous NLS equation and its 
applications in science is meanwhile impossible to summarise. A rather incomplete list of references 
dealing with mathematical aspects of the AL equation is [3,8-46], and [47-69] more generally for 
its hierarchy. Physical applications of the AL equation can be found e.g. in [70-76]. In this work, 
we consider more generally matrix versions of the (continuous and discrete) NLS equations and 
present corresponding solutions. Matrix generalisations, including vector versions, of the scalar 
continuous NLS equation have been studied in particular in [3,77-101]. A brief summary of the 
physical relevance can be found in [3]. In particular, a (symmetric) 2x2 matrix NLS equation 
turned out to be of relevance for the description of a special Bose-Einstein condensate (with atoms 
in a spin 1 state) [102-115]. Semi-discrete matrix NLS equations appeared in [3,116-129], a full 
discretisation has been elaborated in [130] and a dispersionless limit studied in [127]. 

In section 2 we introduce the general framework. Section 3 then presents bidifferential calculi 
for the NLS system and its discrete versions. Section 4 derives exact solutions, by application 
of the method of section 2. The solutions are parametrised by matrix data that have to solve a 
Sylvester equation.^ In section 5, these results are translated into a decomposed form, where the 
NLS systems attain a more familiar form. 

Section 6 deals with a complex conjugation reduction which applies to the continuous as well as 
the discrete NLS systems. Quite surprisingly, this has apparently not been studied previously for 
matrix NLS equations. Our analysis in particular addresses regularity and asymptotic behaviour 
of matrix soliton solutions. Section 7 treats the more familiar Hermitian conjugation reduction 
which, however, only works for the continuous NLS system, at least without severe restrictions or 
introduction of non- locality (cf. appendix D). Section 8 contains some concluding remarks. 

Some supplementary material has been shifted to appendices. Appendix A briefly treats the 
considerably simpler example of the matrix KdV equation, in particular to provide the reader with 
a quick access to the methods used in this work, but also in order to stress the similarities even 
in details of the calculation. Appendix B establishes a relation between the semi-discrete matrix 
NLS (AL) system and the matrix modified KdV equation. Appendix C demonstrates that our 
fully discrete matrix NLS system reduces in the scalar case to the corresponding equation in [5]. 
Appendix D presents a non-local reduction of the semi-discrete NLS system. Appendix E recovers 
Lax pairs for the three NLS systems, starting with the respective bidifferential calculus. Finally, 
Appendix F contains some remarks on the continuum limit of the fully discrete matrix NLS system 
and the solutions obtained in this work. 



See e.g. [131-133] for results on the Sylvester equation. A vast literature deals with its solutions and applications. 
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2 The framework 



Definition 2.1. A (complex) graded algebra is an algebra^ Q, that has a direct sum decomposition 

J7 = 0Q^ (2.1) 

T->0 

into complex vector spaces Q^' , with the property 

n'' c n''+' . (2.2) 

Furthermore, A := 0.^ is a subalgebra and, for r > 0, il^ is an ^-bimodule. 

Definition 2.2. A bidifferential graded algebra or bidifferential calculus is a graded algebra 

equipped with two graded derivations d, d : — )• of degree one (hence dCl^ C 0^+^, dO^ C r2''+^), 
with the properties 

dod = 0, dod = 0, dod + dod = 0, (2.3) 
and the graded Leibniz rule 

d(x X') = (dx) X' + (-l)^' X dx' , d(x xO = (dx) X + ("l)^' X dx' , (2.4) 
for all X G ^' and x' £ 

For an algebra A, a corresponding graded algebra is given by 

n = A®f\{C'^), (2.5) 

where /\{C'^) denotes the exterior algebra of C^.^ In this case it is sufficient to define the maps d, d 
on A. They extend in an obvious way to such that the Leibniz rule holds. Elements of fXiC^) are 
treated as "constants". This structure underlies all our examples. In this work, C1X2 will denote 
a basis of f\\C^). 

Given a unital algebra B and a bidifferential graded algebra (0,d, d) with = B, the actions 
of d and d extend componentwise to matrices over B. A corresponding example is provided in 
Appendix A, the matrix KdV equation. By application of results stated below we recover in a 
strikingly simple way a class of KdV solutions derived previously via inverse scattering theory (see 
in particular [134]). A treatment of NLS-type equations requires a generalisation, however. Here d 
and d involve a structure that depends on the size of the matrix on which they act. 

The space of all matrices over B, with size greater or equal to that of uq x no matrices,^ 

Mat„o(i3) := Mat{n',n,B), (2.6) 

n' ,n>no 

has the structure of a complex algebra with the usual matrix product extended trivially by setting 
AB = whenever the dimensions do not match. Extending A = Ma,tna{B) to a graded algebra (2.5), 
we shall require that d and d respect the size of matrices in the sense that d(Mat(n', n, B)(8)/\(C^)) C 
Mat{n' ,n,B) (8) /\(C^)), and correspondingly for d. The following result originated in [135] and 
appeared in [1] in a more general form. 

^All algebras in this work will be (assumed to be) associative. 
^More generally, one may consider the exterior algebra of C*^, fc > 1. 

*In some cases it may turn out to be necessary to further restrict the values of n and n'. In typical examples we 
can set no = 1, but the NLS case requires no = 2 (see section 3). 



4 



Theorem 2.3. Let (r2,d,d) be a bidifferential graded algebra witit' Q = A® A(*C^) ^'^'^ — 
Matno(^), f'^'^ some no S N. For fixed n,n' > tiq, let X G Mat(n, n,;S) and Y e Mat(n',n,;S) &e 
solutions of 

dX = (dX) P , dY = (dY) P , (2.7) 

andP 

RX-XP = -QY, (2.8) 

with d- and d-constant matrices P,R G Mat(n,n, ;B), and 

Q = VQU, (2.9) 

where U S Mat(m', ;B), V E Mat(n,m,;B) and Q £ Mat{m,m' ,B) are also constants of the 
calculus. IfX. is invertible, then 

(j) := UYX^^V G Mat(?n,', m, B) (2.10) 

satisfies 

d0 = (d(^) Q + d^ with = UYX-1rV. (2.11) 
Furthermore, (by application of d) (j) solves 

dd(t> = d(t)Qd(f> . (2.12) 
Proof: As a consequence of (2.7) and (2.8), we find (cf. [1]) that 

* := YX"^ G Mat(n',n,B) 

solves 

d* = (d*)Q* + d(*R) . 

Multiplying from the left by U and from the right by V, resolving Q according to (2.9), and using 
the graded Leibniz rule for d and d, we obtain (2.11). □ 

The theorem reveals an integrability feature of the nonlinear equation (2.12). From solutions of 
a system of n x n, respectively n' x n, matrix linear equations, it generates a large class of solutions 
of the m' X m matrix equation (2.12). Since n,n' can be arbitrarily large, in this way we obtain 
infinite families of solutions of (2.14), where the members depend on a number of parameters that 
increases with n and n'. In this work we will restrict our considerations to the case where 

m' = m , Q = Im ■ (2.13) 

Accordingly, our central equation will be 

dd(p = d<j) d(j) where </> G Mat(m, m, 5) . (2-14) 



^More general choices of Q, are actually allowed. 

®(2.8) has a form analogous to that of a Sylvester equation. The ingredients of the latter are just matrices over 
C, however. Nevertheless, we shall see that in our examples (2.8) actually results in an ordinary Sylvester equation. 
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Remark 2.4. Concerning the general case, if Q admits a factorisation 

Q = a( ^r>,(m'-r) \ ^ 

\ ^(m—r)xr '^{m—r)x{m'~r) / 

into d- and d-constant matrices, where A and B are invertible, then (2.12) is equivalent to a 
(inhomogeneous) linear extension (cf. [136]) of (2.14) (where m has to be replaced by r), after the 
redefinition B(j)A i— t- (/). This holds in particular if Q is a constant complex matrix (Smith form 
decomposition, see e.g. [137]), and if all such matrices are d- and d-constant. Our methods can 
also be applied to such a linear extension. In case of the bidifferential calculi in section 3, d- and 
d-constancy of a matrix places additional restrictions on a constant matrix, so that (2.12) is not 
in general equivalent to a linear extension of (2.14). All this opens additional possibilities that 
will be disregarded in this work, however. A matrix Q different from an identity matrix effects a 
deformation of the ordinary matrix product, which indeed offers interesting applications [138]. □ 

Many (e.g. in the sense of the inverse scattering method) integrable partial differential or 
difference equations can be recovered from the universal equation (2.14), by choosing a suitable 
bidifferential graded algebra. At first sight it seems that (2.14) can only produce equations with 
quadratic nonlinearity. But (p need not directly correspond to the independent variable of the 
equation of final interest. Indeed, the main examples in this work, the NLS systems, possess a 
cubic nonlinearity. The (anti-) self-dual Yang-Mills equation (more precisely, a gauge reduced 
version of it), which is well-known as a source of many integrable partial differential equations, is 
a special case of (2.14) [1]. Whereas the KP equation and its hierarchy do not quite fit into the 
self-dual Yang-Mills framework, they do fit well into our generalised framework. The generalisation 
(2.14) covers much more, however, and in particular integrable partial difference equations, i.e. 
discretised versions of continuous integrable equations. 

In all examples that will be considered in this work, the general solution of the linear equations 
(2.7) can be expressed in the form 

X = Ao + AiH, Y = Bo + BiH, (2.15) 

where the matrices Ao,Ai,Bo,Bi are d- and d-constant, and H is not. Assuming that Aq is 
invertible, we can rewrite the solution (j) determined by the theorem as 

= U (Bo + Bi H)(I - KH)~^V = U3 (I - KH)-^V + UBqV (2.16) 

in terms of the d- and d-constant matrices 

K = -A^^Ai, U = U(Bi + BoK), V = A^ . (2.17) 

The last summand in (2.16) is d- and d-constant. Since such a term can always be added to any 
solution of (2.14), we shall neglect it. Next we insert (2.15) in (2.8) and separate the d- and d- 
constant summands from those that depend on H. Assuming that both parts have to be satisfied 
separately (which is true in our examples), the first part becomes R = (AqP — VUBq) Aq"^. 
Eliminating R from the second, then yields 

[P,K3]=VUS. (2.18) 

We summarise our results. 
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Corollary 2.5. Let {^[l = /\(C^),d, d) be a bidifferential graded algebra with A = MainQ^B) 
for some hq € N. Let H G M.at{n,n,B) be any solution of 

d3 = (dH)P (2.19) 

and (2.18), with d- and d-constant matrices P,K € Mat{n,n,B), V G Mat{n,m, B) and U G 
Mat(m,n,i3). Then 

(j) = JJ S {I - KSy^V (2.20) 

solves (2. 14-), provided that I — K3 is invertible. □ 

Having determined S (which should not be d- and d-constant) in a concrete example, the 
remaining task is thus to find, for given d- and d-constant matrices U and V, a d- and d-constant 
matrix K such that (2.18) holds. 

Remark 2.6. The construction of exact solutions according to theorem 2.3 hardly exhausts the 
set of all solutions of (2.14) (which is (2.12) with Q = /„). In examples like matrix KdV (see 
appendix A) and the matrix NLS systems treated in the main part of this work, the latter equation 
is equivalent to di^ = (d(/>) 4> + dO with an arbitrary matrix function 9. But (2.11) chooses from 
a certain family of matrix functions, which clearly means a restriction to a subclass of solutions. 
On the other hand, it should be noted that the expression for 6 in (2.11) involves matrices of 
arbitrarily large size and, with suitable technical assumptions, we may allow them to be infinitely 
large, corresponding to operators on a Banach space (see also [139-142]). It is by far not evident 
what kind of solutions can be reached in this way. □ 

Remark 2.7. Theorem 2.3 can be generalised by replacing (2.7) with the weaker equations 

dX = (dX)P + X*, dY = (dY)P + Y*, (2.21) 

with any G $7^ of size n x n, and dropping the requirement that the n x n matrix P is d- and 
d-constant. P and ^ are constrained as a consequence of the above system, however. Indeed, 
acting with d on (2.8), using the Leibniz rule, (2.21) and (2.8) again, we obtain 

dP = (dP)P + P*-*P . (2.22) 

Acting with d on (2.21), using the Leibniz rule and dd = — dd, again (2.21) and then (2.22), leads 
to 

d* = (d*)P-*^. (2.23) 

Setting 

* = (dG - (dG) P) , (2.24) 

with an inver tible nxn matrix function G, then = d^G = d(G* + (dG)P) = G (d* - (d*) P + 
*2), by use of (2.22), shows that (2.23) is satisfied. In terms of P' = GPG~^ (2.22) then reads 

dP' = (dP')P'. (2.25) 

If for any solution (P, ^) of (2.22) and (2.23) the equation (2.24) has a solution G, it follows that 
^ can be set to zero without restriction of generality. A restriction seems to remain nevertheless 
in theorem 2.3, since there we solved the non-linear equation (2.25) trivially by requiring that P' 
(renamed to P) is d- and d-constant. In case of other solutions of (2.25), the linear equations (2.7) 
involve variable coefficients and also the constraint (2.8) becomes considerably more complicated. 
We will therefore not consider this generalisation further in this work. □ 
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3 NLS systems 



3.1 A class of bidifferential calculi 

Let A = Mat2(.S) and O = A® /\{C'^). For every n > 2 we choose a matrix Jn with the properties 

4 = In, Jn^±In, (3.1) 

where In denotes the n x n identity matrix. Then 

e2(/):=^[J,/]:=^(Jn'/-/Jn) for / G Mat(n', ?i, (3.2) 

defines a derivation of A. Setting 

df = eiif) Ci + e2(/) C2 , df = ei(/) Ci + e2(/) C2 , (3.3) 

then determines a bidifferential calculus on 0, if ei,ei,e2 are commuting derivations of A and 
also commute with 62- In the following, 61,61,62 will be related to partial derivative or partial 
difference operators. In order to be d- and d-constant, a matrix must then not depend on the 
respective variables. We shall see that, with appropriate choices, (2.14) reproduces a continuous, 
semi-discrete, respectively fully discrete NLS system. 

3.2 Continuous NLS 

Let B be the space of smooth complex functions on M^. In (3.3) we choose^ 

ei = dx, 61 = -i dt, €2 = dx, (3.4) 
where dx is the partial derivative with respect to x, and i = (2.14) takes the form 

-^[J,(t't] = (pxx + ^[(l)x,[JA]], (3.5) 

where J = Jm- Writing 

cf) = J if , (3.6) 

and decomposing the m x m matrix ip as follows into the sum of an "even" and an "odd" part, 

(p = V + U, V = ^{ip + J^J), U = ^{ip - JifJ) , (3.7) 

we have 

JV = VJ, JU = -UJ . (3.8) 

Hence (3.5) splits into 

UUt+Uxx+'DxU+UVx = (3.9) 
and Vxx = —{U'^)x- The latter can be integrated to 

Vx = -U'^ + C, (3.10) 



^Despite of our notation, in some important applications (notably in optics) x plays the role of a time coordinate 
and t that of a physical space coordinate. 
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where C is an even matrix independent of x. This can be used to eliminate D in (3.9). We conclude 
that, for the chosen bidifferential calculus, (2.14) is equivalent to 

lJUt+Uccx-2U^ +CK +KC = 0, Vx = -U"^ + C . (3.11) 

Without restriction of generality, we can set^ C = 0, so that 

lJKt+Uxx-2U^ = 0, Vx = -U'^. (3.12) 

This will be called (m x m matrix) continuous NLS system. 
The C2-part of (2.11) (with Q = /„,) is 

= ^[J,(^] J^+^J[J,e], (3.13) 

which implies Vx = —U^- As a consequence, the solutions generated via theorem 2.3 actually have 
vanishing C.^ 

3.3 Semi-discrete NLS 

Let Bq be the algebra of functions on Z x M, smooth in the second variable. We extend it to an 
algebra B by adjoining the shift operator S with respect to the discrete variable (the discretised 
coordinate x). In (3.3) we choose 

ei(/) = [§,/], ei(/) = -i/, e2(/) = -[S-i,/], (3.14) 

where / = df /dt. In the following we use the notation = SfS~^, = E^^fS. Setting 

= JifS^^ , (3.15) 

(2.14) becomes 

^[J,ip] + ip+ -2^ + ip- + - ^)J[J, 99] - i[J, ifjJi^ - if-) = , (3.16) 

where 93 can be restricted to be an element of Mat(m, m, Bq). With the decomposition ip = V + 14 
as used above, we find 

V+ - V = -U+U + C , (3.17) 

where C is even and does not depend on the discrete variable, and 

iJU+U+ -2U+U- - U^U"^ -U^U- + {C,IY} = . (3.18) 

With the above choice of bidifferential calculus, (2.14) is thus equivalent to (3.17) and (3.18). Again, 
C can be eliminated by redefinitions, and we have 

iJU + -2U +U- -U+U^ -U'^U- = 0, V+-V = -U+U. (3.19) 

These two equations constitute a semi-discretisation of (3.12) and will be called (m x m matrix) 
semi-discrete NLS system. By inspection of (2.11), the solutions generated via theorem 2.3 have 
vanishing C. 



^For U' = TUT'^, where Tt = -iJCT, and V = V ~ Cx, (3.12) holds without the terms mvolving C. The 
modification of the NLS equations in sections 6 and 7, caused by a nonvanishing C, in the scalar case reproduces a 
"background term" which turned out to be convenient in order to study finite density type solutions (see e.g. (1.5) 
in [2]). 

^We actually lost the freedom of the choice of an even matrix C in the step from (2.11) to (2.12), which would 
still allow the addition of a d-closed 1-form that is not d-exact, namely JC (2 in the case under consideration. 
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3.4 Fully discrete NLS 



Here we make the following choice in (3.3), 

ei(/) = [TS, /] , ei(/) = -i [T, /] , e2(/) = f] , (3.20) 

where T is the shift operator in discrete "time". Here B is the algebra of functions Bo on Z x Z, 
extended by adjoining S and T. In the following, we also set /+ = T/T~^ and /_ = T~^/T. 
Writing 

(p = J^S-^ , (3.21) 

(2.14) reads 

^[J,(p+ + (V5+ - f)+ -{f- + \{^X - ^)J[J,ip+\ - \[J,'P\J{^+ - V'") = 0, (3.22) 

and if can be restricted to be an element of Mat(m, m, So). With the decomposition if = D + U , 
we obtain 

iJ{U+-U) + {U+ - U)+ -{U- U~) -U+Ul- U^U- + {V+ - V)U++U{V+-V) = 0, 
V+ -V = -U+U , (3.23) 

which will be called {m x m matrix) fully discrete NLS system. We again dropped terms involving 
an even matrix C that does not depend on the discrete variable x. Again, inspection of (2.11) shows 
that solutions generated via theorem 2.3 have in fact vanishing C and are therefore solutions of 
(3.23). 



3.5 Normal form of the matrix J 

By an application of the Jordan decomposition theorem, J is related by a similarity transformation 
to a diagonal matrix with diagonal entries ±1. Since we can apply this similarity transformation 
to the above equations, without restriction of generality we can choose 

with any choice mi,m2 G N such that mi + m2 = m. Note that, for m > 2, the bidifferential 
calculus depends on the splitting of m into two summands. This means that there are actually 
different bidifferential calculi on Mat(m, m, B), and different choices correspond to different NLS 
systems. 



3.6 Complex conjugation reduction 

Let J satisfy 

J* = J, (3.25) 

where * denotes complex conjugation, and let F be an invertible constant m x m matrix with the 
properties 

JT = -TJ, r* = r-^. (3.26) 
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Then the transformation 



U^er~^U*r with e = ±l, V ^ r^^V*r , (3.27) 

leaves ah three NLS systems invariant, including (3.8) (which applies to all three cases). Hence a 
consistent reduction is given by 

w = erur-^, v* = rvT-^. (3.28) 

Transforming J to its normal form (3.24), the first of (3.26) restricts F to off-block-diagonal form. 
The second of conditions (3.26) then requires that m is even and 

mi = 7712 ■ (3.29) 

We will actually choose 




(3.30) 



3.7 Hermitian conjugation reduction 

The continuous NLS system with 

/ = J, (3.31) 

where ^ denotes Hermitian conjugation (i.e. transposition combined with complex conjugation), is 
invariant under 

U^eU^ with e = ±l, V ^ V'^ . (3.32) 

Hence 

= eU, = V (3.33) 

is a consistent reduction. This does not work as well for the matrix discrete NLS systems, since 
^ is an anti-involution (whereas * is an involution). Applying it to the nonlinear term in the first 
of equations (3.19), we have {U+U^ +U^U~)^ = + U^-U^^. Clearly the matrices in this 

expression will not commute in general. In order to achieve that the equations obtained from (3.19) 
by Hermitian conjugation become equivalent to the original equations, we can either 

(1) require that commutes with lA^, while keeping (3.33), or 

(2) replace (3.33) hyU{x,ty = eU{-x,t) and P(x,t)'f = -V{-x,t) (see also [116]). 

We will not discuss the first case further in this work, but only refer to [3,92, 118, 123, 124, 130], see 
also remark 7.5. The second case leads (in a decomposed form in the sense of section 5) in general 
to a non-local equation, see appendix D. 

4 Solutions of the matrix NLS systems 

A central result of this section is expressed in the next proposition. In the following, I = /„ denotes 
the 71 X 77 identity matrix. 
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Proposition 4.1. 

= U(3-^ - KHK)-^ V, P = UHK(H-^ - KSK)-i V (4.1) 

solves^^ the continuous, semi- discrete, respectively fully discrete mxm matrix NLS system (involv- 
ing the m X m matrix J with = Im) if 





J / ±1 


S 


U 


V 


K 


size 


n X n 


n X n 


m X n 


n X m 


n X n 



are constant matrices such that 

j2 = i^ JS = SJ, JU = -UJ, JV = VJ, JK = -KJ, (4.2) 

and 



NLS 


Sylvester equation 




continuous 


SK + KS = VU 




semi- discrete 


S-^K - KS = VU 


gx g-it (S+S-1-2I) J 


fully discrete 


S~^K - KS = VU 


S- [(I + i J (I - s-i))(i - i J (I - s))-i]* 



where x £ in the semi- and fully discrete NLS case, t £ Z in the fully discrete NLS case, and 
x,t otherwise. Of course, we need S invertible in the discrete cases. □ 

The matrices S,U,V,K are subject to a Sylvester equation. Given matrix data (J,S,U,V), 
the problem is to find a solution K of the Sylvester equation. Let o"(S) denote the spectrum of S, 
i.e. the set of its eigenvalues. A sufficient condition for the existence of a solution of the respective 
Sylvester equation is 

(t(S) n o"(— S) = continuous NLS , , 

cT(S)nc7(S-i) = discrete NLS ^ ^ 

(see e.g. theorem 4.4.6 in [133]). Moreover, the solution is then unique, and as a consequence K 
automatically satisfies JK = — KJ. 

There is a certain redundancy in the matrix data, since different matrix data can determine the 
same solution, but we will narrow this down to a considerable extent (section 4.2, see also [143]). 

We note that the NLS equations can determine P only up to addition of an arbitrary m x m 
matrix T>q that does not depend on x (and Pq has to commute with J). We shall therefore 
understand the expression for T> in (4.1) always modulo addition of an arbitrary constant matrix 
commuting with J. To put it another way, we simply drop any constant summand that arises in 
an expression for D, and we actually did this in order to arrive at the expression in (4.1). 

In section 4.1 we derive the results formulated in proposition 4.1. Some symmetries of the 
solution formulae in proposition 4.1 are revealed in subsection 4.2. Subsection 4.3 formulates a 
superposition rule on the level of matrix data. Furthermore, section 4.4 specialises to the complex 
conjugation, respectively Hermitian conjugation reductions. 

^°More precisely, here we assume that the inverses appearing in these formulae exist. This will also be done in the 
formulation of results that originate from this proposition. A considerable part of this work deals with this regularity 
problem, however. 
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Remark 4.2. The results formulated in proposition 4.1 will be derived below by application of 
corollary 2.5. The latter actually allows to replace H in proposition 4.1 by CS with any constant 
matrix C that commutes with J and S. But the line of arguments that led us from theorem 2.3 
to corollary 2.4 in fact shows that such a matrix can be eliminated by redefinitions of the matrix 
variables entering the solution formula. Conversely, in special cases (see e.g. section 6.2) it will be 
possible and convenient to shift all the freedom in the choice of U and V into such a matrix C. In 
other cases it will be convenient to move at least part of the freedom from U and V to a matrix 
C. Although the corresponding reparametrisation transformation (see section 4.2.3) requires some 
restrictions that result in an invertible C (commuting with S), the aforementioned ensures us that 
this restriction for C is not really required in order to obtain an NLS solution. □ 

4.1 Proof of proposition 4.1 

In the following we start with the bidifferential calculi on 17 = Mat2(;S) C?) /\{C^) associated in 
section 3 with the continuous NLS, semi-discrete NLS, respectively fully discrete NLS system. 
Then we apply theorem 2.3 in the form specified in corollary 2.5. For fixed n we write J for J„. 
The requirement that the matrices P, U, V have to be d- and d-constant now means that they 
must not depend on the independent (continuous, respectively discrete) variables, and they have 
to satisfy 

[J,P]=0, JU = UJ, JV = VJ. (4.4) 

Any n X n matrix A decomposes into the sum of an "even" and an "odd" part, 

A = Ae + Ao where JAg = AgJ , JAq = -AqJ . (4.5) 

Such a matrix is d- and d-constant iff it is even and constant (i.e. independent of x and t). A 
constant n x n matrix 5] that satisfies 

I]2 = I, JS = -S;j, (4.6) 

can be used to convert odd into even matrices and vice versa. 

Continuous NLS. We write 

P = JS, (4.7) 
with a new object S. Then dX = (dX) P = (dX) J S becomes 

-iXt = X,SJ, X, = i[J,X]JS. (4.8) 

The general solution is 

X = Ao + AiS3 with g ^ ^-xS-its^j ^ ^^^ g^ 

Here Aq, Ai are constant even matrices (hence d- and d-constant). Assuming that Aq is invertible, 
application of theorem 2.3 reduces to an application of corollary 2.5. In terms of the redefined 
matrices 

K = K5:, U = US, (4.10) 
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and with H = S3, (2.18) takes the form 

SK + KS = VU, (4.11) 

Note that JV = VJ, but JU = -UJ. (2.20) and (3.6) lead to 

99 = U3 (I - KH)-^V = UH (I + KH)(I - {KSfy^V . (4.12) 

The last expression is convenient in order to split ip into an even and an odd part. In this way we 
obtain (4.1). 

Semi-discrete NLS. Writing 

X = SX, P = JS~^P, (4.13) 
we find that dX = (dX) P = (dX) J S^^P is equivalent to 

-iX= (X+-X)JP, X-X~ = ^[J,X]JP, (4.14) 
with the general solution 

X = Ao + Ai5]H, (4.15) 

where 

g^gxg-ij^(s)j^ a;(S) = S + S-^ -21, S = (I + P)-^ (4.16) 

and X G Z. Using the same redefinitions (4.10) as in the continuous NLS case, (2.18) becomes 

S"^K-KS = VU, (4.17) 
and (4.1) is obtained in precisely the same way as in the continuous NLS case. 

Fully discrete NLS. Using again (4.13), we find that dX = (dX) P is equivalent to 

-i(X+ -X) = (X+ -X) JP, X- X- = ^[J,X] JP . (4.18) 

Assuming that I — i JP is invertible, it follows that the even part of X is constant. If furthermore 
S = (I + P)~-^ and [I — i J(I — S)]~^ exist, the above system yields the following equations for the 
odd part of X, 

X+ = XoS, (Xo)+ = Xo[I + iJ(I-S-i)][I-iJ(I-S)]-^ (4.19) 
Hence we obtain the solution 

X = Ao + Ai E H , (4.20) 

with constant even coefficient matrices and 

H = S^'[(I + iJ(I-S-i)][I-iJ(I-S))-Y x,t£Z. (4.21) 

Again, assuming Aq invertible and adopting the definitions (4.10), we arrive at (4.1) and (4.17). 
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4.2 Transformations of the matrix data determining solutions 

4.2.1 Similarity transformations 

A transformation 

J ^ MJM^^ , K ^ MKM~^ , S ^ MSM~^ , U ^ UM"^ , V ^ MV , (4.22) 

where M is any constant invertible n x n matrix, obviously leaves l/( and D in proposition 4.1 
invariant, and also the corresponding Sylvester equation. It can thus be used to achieve that J 
takes the form 

'o -I)' 

with ni,n2 G N such that ni + n2 = n. The remaining similarity transformations such that 

[J,M] = (4.24) 

can then be exploited to diminish the redundancy in the "parametrisation" of the solutions by the 
set of matrices S,U, V. In particular, it can be used to transform S to Jordan normal form. 

4.2.2 Reflection symmetries 

If K is invertible and 11 a constant n x n matrix with 

jn = -nj , (4.25) 

then the transformation 

Sh^s', Uh^-UK-^n-^ Vh^nK-^v, Kt^nK^^n-^ (4.26) 

with 

_ J — nSn~^ continuous NLS , . 

~ \ nS'^n-i semi- and fully discrete NLS ^ ' ^ 

implies 

Hh^ns-^n^^ (4.28) 

and leaves V( in proposition 4.1 and also the corresponding Sylvester equation invariant. We refer 
to such a transformation as a reflection. Since 

V^V (HK)~i (H-i - KHK)-^ V = P + U K"^ V , (4.29) 

T> only changes by a constant matrix, which can be dropped as pointed out above. Hence a 
reflection is a symmetry of the solutions given in proposition 4.1 for all the NLS systems. A 
symmetry transformation of this kind for continuous NLS solutions already appeared in [143]. 

4.2.3 Reparametrisation transformations 

Let A,B be any constant invertible nxn matrices that commute with J and S, i.e. 

[J, A] = [J, B] = [S, A] = [S, B] = . (4.30) 

Then 

Sh^S, Uh^UB"\ Vh^A^^V, Kh^A"^KB"\ Hh^BAH, (4.31) 

is a reparametrisation of the solution (4.1). This is helpful in order to reduce the number of 
parameters of U and V (see in particular section 6.2). 
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4.3 Superposition of solutions 

Let (Jj, Sj, Uj, Vj), i = 1,2, be data that determine two solutions {Ui^Vi) of one of the three 
systems (3.12), (3.19) or (3.23). Let Kj, i = 1,2, be corresponding solutions of the Sylvester 
equations. Then 

determine new solutions if K12 and K21 satisfy 

J1K12 = -K12J2 , J2K21 = -K21J1 , (4.33) 

and solve the respective equations below. 
Continuous NLS: 

S1K12 + K12S2 = V1U2, S2K21 + K21S1 = V2U1 . (4.34) 

These equations possess unique solutions if cr(Si) n (7(— S2) = 0- 
Semi-discrete and fully discrete NLS: 

Sr^Ki2-Ki2S2 = V1U2, S^iK2i - K21S1 = V2U1 . (4.35) 

These equations have unique solutions if (T(Si)n(T(S^^) = 0. (4.3) applied to S = block-diag(Si, S2) 
contains the abovementioned spectrum conditions and is therefore a sufficient (but not necessary) 
condition for the above superposition principle to generate again a solution. Let us list some fairly 
obvious properties of superpositions. 

1. The superposition of two solutions is commutative. Note that 

is a similarity transformation that exchanges the matrix blocks corresponding to the two 
solutions in the superposed data. 

2. Superpositions are associative. 

3. A similarity transformation of the matrix data corresponding to a constituent solution extends 
to a similarity transformation of the superposition, hence remains an equivalence transfor- 
mation. This is simply based on the fact that a similarity transformation of the superposed 
data, with an invertible matrix of the form 

covers arbitrary similarity transformations of the constituent data. 

4. A reparametrisation transformation of a constituent extends to a reparametrisation of the 
superposition, hence remains an equivalence transformation. (Choose A, B block-diagonal in 
(4.30) and (4.31).) 
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Recall that, for any data (J, S, U, V) that determine a solution, via a similarity transformation 
J can be brought to standard form (4.23) and S simultaneously to Jordan normal form. We call 
the data (J,S,U, V) simple if at least one of the blocks of S, corresponding to a block of J, is a 
single Jordan block. For any solution given by data (J,S,U, V) that are not simple, we have the 
block structure in (4.32). These data have to solve the Sylvester equation, hence in particular there 
are matrices Ki,K2 solving its diagonal blocks. It follows that (J,S,U, V) decomposes into data 
(Jj, Sj, Uj, Vj), i = 1,2, that determine solutions. As a consequence, for any solution obtained 
from non-simple data, the simple subdata also determine solutions. In this sense, the solution can 
always be regarded as a superposition of solutions. 

Though reflection symmetries are not quite compatible with the structure of a superposition, 
we still have the following result (based on arguments in [143]). 

Proposition 4.3. Let (4-3) hold for S = block-diag(Si, S2). A reflection symmetry of (only) one 
of the two matrix data sets can then be extended in such a way that it is also a symmetry of the 
superposition solution. 

Proof: We consider a reflection symmetry on the second data set. Hence we have to assume that 
K2 is invertible. Writing Zj = H-^^ — KjHjKj, so that Ui = UjZ~^Vj, we have U = UZ~^V with 



Z 



Zi — K12H2K21 — K1H1K12 — K12H2K2 
-K2H2K21 — K2i3iKi Z2 — K21H1K12 



According to the reflection transformation rules (4.26), we have 

_ Si 
^ " I S'2 

where S2 is defined for the three NLS cases in (4.27). The other transformed quantities U2, V2,K2 
are given by (4.26). Corresponding transformation rules for K12 and K21 should then be read off 
from (4.34), respectively (4.35). It turns out that also Ui and Vi have to be transformed. In fact, 
setting 

„. . . V' = B V „th A = ( _ ) , B = ( J -K-K^ 

and 

^, _ / Ki - Ki2K^^K2i -Ki2K^^n-i \ _f Ki K12 \ . 

^ ~{ nK2iK2i nKs^n-i J ~ \ o -n J ^' 

one obtains Z' = BZA and hence U' = U' Z'~^ V' = U. Furthermore, we find that 

i.e. T> only changes by an irrelevant constant summand. □ 

If we superpose the same data, we obtain nothing new. More generally, we have the following 
result. 



"^One can still refine this proposition in order to weaken the assumptions. 
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Proposition 4.4. Let (Jj, Sj, Uj, Vj), i = 1,2, be two data sets that determine solutions of an NLS 
system according to proposition 4-1 <ind let (4-3) hold for S = block-diag(Si, 82). Suppose there 
are matrices A,B such that U2 = UiB, V2 = AV^, ASi = S2A and SiB = BS2. Let further 
C = I + BA be invertible. Then the superposed data (4-32) determine a solution that can also be 
obtained from the data (Ji, Si, UiC, Vi) (i.e. simply by a redefinition of \Ji). 

Proof: We consider the continuous NLS case only. The discrete cases are treated in a similar way. 
The assumptions together with the respective Sylvester equations for the data (Jj, Sj, Uj, Vj) imply 

S2 (K2 - A Ki B) + (K2 - A Ki B) S2 = . 

Since the assumption (4.3) in particular implies (7(82) H (7(— S2) = 0, the last equation implies 
K2 = AKi B. For (4.32) to determine a solution, (4.34) has to be satisfied, and thus 

Si (K12 - Ki B) + (K12 - Ki B) S2 = , S2 (K21 - AKi) + (K21 - AKi) Si = . 

Since (4.3) also implies (t(Si) n(T(— S2) = 0, these equations have the unique solutions K12 = Ki B 
and K21 = AKi, hence 

K = ( ^ 1 Ki f I B 



A 



This in turn implies 



W:=I-HKHK = I- I ^ HiKiCHiKi ( I B), 

and thus, assuming that C = I + BA is invertible and introducing Hi = C Hi, 

W-i = I+({)c-i(Wr^-I)(l B), 

where Wi = I — Hi Ki Hi Ki. With straightforward computations (4.1) now leads to 

U = Ui (Hr^-KiHiKi)"'vi, P = UiHiKi (h^^ - Ki Hi Ki) Vi , 
from which our assertion is easily deduced. □ 

4.4 Imposing reduction conditions on the solutions 
4.4.1 Complex conjugation reduction 

Let r be an n X n matrix with the properties 

jr = -rj , r* = r^^ . (4.38) 

Choosing the normal form (4.23) of J, this implies that T is off-block-diagonal. The last condition 
in (4.38) requires that n is even and 

ni = n2 . (4.39) 
Imposing the following constraints on the matrices S,U, V, 

s* = rsr-^ , u* = e e' rur"^ , v* = e' rvr-^ , (4.4o) 
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with e' = ±1 and T as in section 3, the Sylvester equations demand 



K* = erKr-\ (4.41) 

and U and T> given by (4.1) indeed satisfy the reduction conditions (3.28). 

The transformations considered in section 4.2 should now be restricted in such a way that the 
reduction conditions are preserved. The matrix M of a similarity transformation (preserving the 
form of J), now also has to satisfy M* = FMr^^. For a reflection symmetry, the matrix 11 has 
to satisfy U* = eTUT-^. In case of a reparametrisation transformation, we have to add the 
conditions A* = TAr-^ and B* = rBr~^ 

4.4.2 Hermitian conjugation reduction 

This only applies to the continuous NLS case. From (4.1) we obtain 

= Vt (St-i - KtstRt)-! , (4.42) 
In order to achieve the reduction condition hi^ = eU with e = ±1, we set 

St = TST"\ U = V^T, (4.43) 
with an invertible n x n matrix T with the properties 

JT = -TJ , = e T . (4.44) 
These conditions again restrict us to the case where 

ni = n2 . (4.45) 

The Sylvester equation now reads 

SK + KS = VV^T , (4.46) 

and is preserved by Hermitian conjugation if 

Kt = eTKT"^. (4.47) 

Then the expression for U in (4.1) indeed satisfies = eU. Furthermore, as a consequence of the 
identity Z~^KH = HKZ~^ where Z = — KHK, the expression for T> in (4.1) satisfies = T>. 

For a similarity transformation with matrix M, which commutes with J, we have to require 
T 1—^ (M^)^-'^TM^. With J given by (4.23) (where ni = 712), this can be used to achieve that 

^=(.1/0 )■ 

After having fixed T, the matrix M of a similarity transformation preserving the Hermitian re- 
duction conditions has to satisfy = TM~^T~^. For a reflection symmetry, 11^ = — eTn~^T~^. 
In case of a reparametrisation, A^ = TBT~^. 
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5 Block-decomposition of the NLS systems and their solutions 

In this section we express the NLS systems in a more famiUar form^^ and translate proposition 4.1 
correspondingly. We choose for J the normal form (3.24) and write 

where p,p,q,q are, respectively, mi xmi, m2 xm2, mi xm2 and m2 xmi matrices. As a consequence 
of (3.8), 

Continuous NLS system. (3.12) becomes the matrix NLS system (see e.g. [3,79,80,144]) 

iqt + qxx-'^qqq = 0, -iqt + qxx-'^qqq = 0, (5.3) 

together with 

Px = -qq: px = -qq- (5.4) 

Semi-discrete NLS system. (3.19) leads to 

p'^ -p = -q^q, P'^ - P= -q^q, (5.5) 

and the matrix Ablowitz-Ladik (AL) system 

iq + q^ -2q + q~ - q+qq- qqq~ = 0, 
—iq + q^ — 2q + q~—q^qq — qqq~ = 0. (5.6) 

Fully discrete NLS system. From (3.23) we obtain 

p+ -p = -q+q, p+-p=-q+q^ (5.7) 

(as in the semi-discrete case) and 

'^{q+-q) + {q'^ -q)+-{q-q~)-{q'^qq)+-qqq~ + {p+-p)q+ + q{p+-p) = o, 
-1(9+ - ^) + (9"^ - ^)+ - (^- ^") - (^^99)+ - ^9^" + (p+ + -p) = 0.(5.8) 

Unlike the situation in the continuous NLS and semi-discrete NLS (AL) cases, here the variables 
p,p cannot be eliminated from the equations for q and q without introduction of non-local terms. 
A reduction of this system appeared in [5] (see also appendix C). Litegrable full discretisations of 
matrix NLS systems appeared in [130]. 



^^See also appendix E for a derivation of Lax pairs for these NLS systems within the bidifferential calculus frame- 



work. 
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Solutions. For J we choose the normal form (4.23). From the properties of the matrices appearing 
in proposition 4.1, we find the following block structures, 

and H is block-diagonal with diagonal blocks H, H. Now proposition 4.1 translates into the following 
statement. 

Proposition 5.1. Let 























U 


U 


V 


V 


S 


S 


K 


K 


size 


mi X 


1712 X ni 


111 X nil 


71,2 X 7712 


111 X Til 


712 X n2 


ni X n2 


71-2 X ni 



he constant complex matrices, with S and S invertible in the discrete NLS cases, and 



NLS 


Sylvester equations 








continuous 


SK + KS = VU 
SK + KS = VU 




^-xS+itS"^ 




semi-discrete 


S-^K -KS = VU 
S-^K - KS = VU 


gx g— itcj(5) 


gx ^tu}{S) 




fully discrete 


S-^K -KS = VU 
S~^K -KS = VU 


5^ VL{Sf 


5^ ^[sy 


X, t G Z 



where 

0J{S) = S + S~^-2In,, 

n{S) = [In,+i{In,-S-^)][In,-i{In,-S)r\ 

n{s) = [/„^-i(/„^-5-i)][/„, + i(4,-5)]-^ 

Then 

q = U {E-^ - KEK)~^ V, q = U {--^ - KEK)~^ V, 
p = UEK (S-i - KEK) y , p = UEK [E'^ - KEK) V , 

solves the respective NLS system. 

5.1 Transformations of the matrix data determining solutions 

The transformations described in section 4.2 decompose correspondingly. 
1. A similarity transformation satisfies (4.24), with J given by (4.23), iff the transformation matrix 
M is block-diagonal, i.e. M = block-diag(M, M). Then we obtain 

S^MSM"^, S^MSM-^, U^UM~^, U^UM"^, 
K^MKM-^, K^MKM-^, V^MV, V^MV. (5.13) 
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(5.10) 
(5.11) 

(5.12) 

□ 



As a consequence, without restriction of generality we can assume that S and S have Jordan normal 
form. 

2. The matrix 11 of a reflection symmetry has to be off-block-diagonal. For later use we consider 
the example 

n^f.? 'i-), (5.14) 



(5.15) 





assuming rii = n2 (and thus n even). The transformation then acts as follows, 
K ^ , U ^ -pUK-\ V ^ p-^K-^V , 



together with 



for continuous NLS, and 



S^-S, S^-S, (5.16) 



S-^ (5.17) 

for semi-discrete and fully discrete NLS. This implies 

E^E~\ Hh^H"^. (5.18) 

3. In case of a reparametrisation transformation, we have A = block-diag(A, A), B = block-diag(-B, B), 
and [S, A] = [S, B] = [S, A] = [S, B] = 0. The transformation is then given by 

S^S, S^S, U^UB'^, U^IJB^^, V^A-'^V, V ^ A'^V 

K ^ A~^KB~^ , k ^ A~^kB~^ , E^BAE, E^BAE. (5.19) 

4. A transformation 

U^SU, tJ^StJ, V^VS-^, V^VS-^, (5.20) 

with invertible constant matrices S,S, has the effect 

q H> SqS~^ , q ^ SqS~^ , p ^ SpS~^ , p ^ SpS~^ . (5.21) 

In contrast to the previous transformations, this is not a symmetry of the solutions. But it may be 
regarded as an equivalence transformation of solutions. 

6 The complex conjugation reduction 

We recall that, for 

mi = 1712, (6.1) 

all three NLS systems admit the complex conjugation reduction introduced in section 3. In the 
following we rename mi to m. Choosing T as in (3.30), the reduction conditions (3.28) take the 
form 

q = eq* with e = ±1 , P = P* , (6-2) 
Here q,q,p,p are all square matrices of size in x m. 
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Continuous NLS. (6.2) reduces the (continuous) NLS system to the following (square) matrix 
NLS equations, 

iqt + qxx-'^eqq*q = 0, Px = -eqq*. (6.3) 

The cases e = — 1 and e = 1 are sometimes referred to as focusing and defocusing, respectively. 
For m = 2 and symmetric q, i.e. q^ = q (where q'^ is the transpose of q), so that q has the form 

qo q~i J 

the matrix NLS equation with e = — 1 describes a special case of a spin-1 Bose-Einstein condensate 
[102-115]. (6.3) is then equivalent to the system 

igi,* + gi,xx + 2(|gi|2 + 2|goP)gi + 2gggli = 0, 
iq-l,t + q-l,xx + 2{\q-i\'^ + 2\qo\^)q^l+2q'^ql = 0, 
i^o.t + 9o,xx + 2(|gip + |go|^ + k-i|^)9o + 2gi g'o 9-1 = 0. (6.5) 

Semi-discrete NLS. In this case (6.2) leads to the (square) matrix AL equations 

iq + q'^ - 2q + q~ - e{q~^q*q + qq*q~) = , p'^ - p = -eq~^q* . (6.6) 

The scalar version of the first equation first appeared in [4-6] and has since been the subject of 
many publications (see the introduction). 

Fully discrete NLS. Now (6.2) yields 

i {q+ -q) + {q^^ - q)+ -{q-q~) {q^q*q)+ - e qq*q~ + {p+ - p)q+ + q{p+ -pT = o, 

p'^ — p = —e q~^q* . (6-7) 

Here p does not decouple from the equation for q. Appendix C shows that, by elimination of p at 
the price of introducing non-local terms in the first equation, one recovers an equation first obtained 
by Ablowitz and Ladik [5] in the scalar case. 

Solutions. We set ni = 712 and rename it to n. Furthermore, we choose 

I. '0 ) • 

Using (5.9), the reduction conditions (4.40) and (4.41) become 

S = S*, U = ee'U*, V = e' V* , K = eK*. (6.9) 
Proposition 5.1 now implies the following result. 

Proposition 6.1. Let S,U,V be constant complex matrices of size n x n, m x n and n x m, 

respectively. In the discrete NLS cases, we assume invertibility of S. Then 

q = ±U {E*~^ - eK*EKyW* , (6.10) 
p = eUE*K* {E-^ - eKE*K*)-^V (6.11) 

solves the mx m matrix continuous NLS equations (6.3), the semi-discrete NLS equations (6.6), 
respectively the fully discrete NLS equations (6.7), if there is a matrix K solving the following 
Sylvester equation and if E is given by the following expression, 
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NLS 


Sylvester equation 




continuous 


SK + KS* = VU 


^-xS-itS^ 


x,t £R 


semi-discrete 


S'^K - KS* = VU 


gx ^—itu}{S) 


xeZ,teM. 


fully discrete 


S~^K - KS* = VU 


5^ n{sY 


x,t G Z 



where 



uj{s) = s + s-^- 2/„ , n{s) = [/„ + i (/„ - s-^)][in -i{in- s)]-^ . 



(6.12) 

□ 



With reference e.g. to theorem 4.4.6 in [133], a sufficient condition for the above Sylvester 
equations to possess a solution K, irrespective of the form of U and V ^ is 



a{S)r\a{-S*) 
a{S*)r\(T{S"^) 



for 



continuous NLS 
discrete NLS 



(6.13) 



(where (j{S) denotes the spectrum of the matrix S). The solution matrix K is then unique. 

Choosing S diagonal, in the scalar focusing case the solutions determined by the last proposition 
contain the familiar solitons. They will be described in a slightly different form in section 6.2. 

Remark 6.2. For given data (S", f7, V) with an n x n matrix S, let i^' be a solution of the Sylvester 
equation in the continuous NLS case. Then it follows that K also solves the Sylvester equation 
with S' = S + i/3/„, /3 G M. Hence there is a new NLS solution obtained by replacing H in (6.10) 
and (6.11) by 



E{x,t-S') = e--^'-i*(^')' = e-'^^^-P'^ E{x - 2/3t,t; S) 
This induces the following Galilean transformation of the NLS solution, 

g'(x,t) = e^^(^-^*)g(x-2/3t,t), 
which is a symmetry of the NLS equation. 



(6.14) 

(6.15) 

□ 



Remark 6.3. If S is symmetric, i.e. S^ = S, and if [/ = 1/^, in the continuous NLS case the 
Sylvester equation becomes the Lyapunov equation^^ 



SK + KS^ = VV^ 



(6.16) 



where K can be chosen Hermitian. Furthermore, the solution q given by (6.10) is symmetric, i.e. 

q^ = q. In particular, this yields solutions of the abovementioned spin-1 Bose-Einstein condensate 

model. For later reference (see remark 6.7), we mention the following result in [148]. 

(6.16) possesses a Hermitian positive definite solution K if and only if the following three conditions 

hold, 

(i) all eigenvalues of S have non-negative real parts, 

(ii) the subspace of C"" spanned by all eigenvectors of S corresponding to eigenvalues with positive 
real part coincides with the controllability subspace associated with the pair (5, V) (which is the 
subspace spanned by the columns of the matrices V, SV, S'^V, . . . , S^~^V), 

(iii) eigenvalues of S with zero real parts are not degenerate (i.e. they are simple roots of the 
minimal polynomial). □ 



^■^More generally, the Lyapunov equation has the form SK + KS^ = Q with a complex matrix Q (see e.g. [132, 
133,145,146]). It plays an important role in control theory [147]. 
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6.1 On the structure of the solutions 



6.1.1 Solving the Sylvester equation 

Without restriction of generality we can assume tliat S has Jordan normal form (since this can be 
achieved by a similarity transformation). The Sylvester equation then has to be satisfied for each 
Jordan block (and in addition we have to solve equations for the off-diagonal blocks of K). Let us 
therefore now consider a matrix S consisting of an r x r Jordan block only, hence 

S = sl + J\f, (6.17) 

where I is the r x r identity matrix and the components of M are given by Mij = Sij-i. Writing 
K = (kij) and 



U = [ Ui • • • Ur ) , V 



(6.18) 



with m-component column vectors Ui,Vi, i = 1, . . . ,r, in the continuous NLS case the Sylvester 
equation takes the form 

{s + s*) kij + ki+ij + kij_i = vjuj i,j = l,...,r, (6.19) 

where kr+ij = ki^ = 0. Assuming s ^ —s* (in which case we know that the Sylvester equation 
has a unique solution), this implies 

^ TUm ' = „ I „* (''^r'^j ~ krj-l) j = 2, . . . ,r , 

hi = ^ ^ {v]ui - fci+i.i) i = r-l,...,l, (6.20) 

which recursively determines the entries kii,kri, i = l,...,r. In fact, (6.19) also determines the 
remaining entries of K recursively. For example, for r = 2 we obtain 

vjui V2U1 v'[u2 V^Ul+V2U2 . r) V2U1 

K = I ~j («+^ ~ is+s")-^ + IJ+Ff 1 211 

I Djltl V2U2 V2UI I ' \ ' / 

s+s* s+s* (s+s*)-^ 

Returning to arbitrary r, we observe that the first column of K vanishes if Ui = 0, and the 
last row of K vanishes if Vr = 0. Inspection of (6.10) and (6.11) shows that in these cases the 
solution is completely determined by the correspondingly reduced matrix data. In the case where 
Ui = 0, these reduced data are obtained by deleting the first row and first column of S, and also 
deleting the first column of U and the first row of V. This also holds in case of the discrete NLS 
versions. We conclude that, for any solution determined by matrix data {S, U, V), where U and V 
are different from zero and 5 is a Jordan block subject to the condition (6.13), there are equivalent 
data (5', U' , V') where S' is a Jordan block (with the same eigenvalue and of size smaller or equal 
to that of S) such that the first column of U' as well as the last row of V are different from zero. 

More generally, for any non- vanishing solution determined by matrix data {S,U,V), with S 
given in Jordan normal form subject to the condition (6.13), we can assume without restriction of 
generality that for any Jordan block of S the first column of the corresponding part of U as well 
as the last row of the corresponding part of V are different from zero. 
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6.1.2 Superpositions 

Given two data sets {Si, Ui, Vi), i = 1, 2, that determine NLS solutions, with corresponding solutions 
Ki of the respective Sylvester equations, we can superpose them, 

where 1^12 and K21 have to solve 

S1K12 + Ki2Sl = ViU2, S2K2i + K2iSl = V2Ui continuous NLS (6.23) 

K12 - S1K12S2 = V1U2 , K21 - S2K21SI = V2U1 discrete NLS . (6.24) 

(6.13) is a sufficient condition for such a solution to exist (which is then unique). Then we have 

/ Zi-e Kl^E2K2i -e KIE1K12 - e Kl^E2K2 

V -e KIE2K21 - e 7^-2*1^1^1 Z2 - e K^^EiKu 

where Zi = E*~^ - eK*EiKi. Its inverse, needed to evaluate q and p more explicitly, can be 
expressed in terms of Schur complements. 

Conversely, given solution data {S,U,V), where S is block-diagonal with at least two blocks, 
then these data can be understood as a superposition (see also section 4.3). 

The following result is simply proposition 4.4 adapted to the case under consideration. It will 
be helpful in the sequel. 

Proposition 6.4. Let {Si,Ui,Vi), i= 1,2, be two data sets that determine solutions of an NLS 
system according to proposition 6.1 and let (6.13) hold for S = block-diag(S'i, S'2). Suppose there 
are matrices A, B such that U2 = UiB, V2 = AVi, ASi = S2A and S\B = B S2. Let further C = 
L + B*A be invertible. Then the superposed data (6.22) determine a solution that can equivalently 
be obtained from the data (Si, C/iC, Vi) (i.e. simply by a redefinition ofUi). 

6.1.3 The use of reflection symmetries 

The following result has its origin in [143], where only the (focusing) continuous NLS equation was 
considered, however. We show that for a solution determined by matrix data via proposition 6.1, 
in general equivalent matrix data (5, C/, V) exist such that all eigenvalues s of S" with Re(s) 7^ in 
the continuous NLS case, respectively \s\ 7^ 1 in the discrete cases, satisfy 

(i) Re(s) > in the continuous NLS case, 

(ii) |s| < 1 in the semi- and fully discrete NLS case. 

In case (i) the matrix S is called positive stable [133]. In case (ii) it is sometimes called stable with 
respect to the unit circle [145]. 

Suppose S'2 is a Jordan block that does not satisfy condition (i), respectively (ii). Writing 

presents the solution data in the form of a superposition, and one has to solve (6.23), respectively 
(6.24), in order to construct the solution. The idea is now to apply a reflection symmetry and to 
use the analogue of proposition 4.3. A reflection symmetry of the form considered in section 4.2.2 
preserves the reduction conditions if 

^ = e/3*. (6.27) 



(6.25) 
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Let us choose /3 = e (so that (3 = 1). The reflection symmetry is then given by 

U2^-eU2K^\ V2^K^^V2, K2^-eU2K^\ (6.28) 

together with 

^ Q> _ j —S2 continuous NLS , . 

^ ^ 2 - I f^uy discrete NLS ' ^ ^ ^ 

This is a reflection in an obvious sense. As in the proof of proposition 4.3, we can construct an 
extension to matrix data that determine the same superposition solution. For the above reflection, 
this is achieved by 



K^l Ki-Ki,_K:^^K2, "'^'i'f^' I . (6.32) 



and 



K:^'K2i ei^2 

The only assumption needed for this to work is that the Sylvester equation with data (5*2, C/2) V2) 
has an invertible solution K2- 

In case of a Jordan block, the result of the reflection (6.29) has no longer the standard Jordan 
form. But with the help of a similarity transformation one can restore the usual Jordan structure 
(without changing the eigenvalues). 



6.2 Solutions of the scalar continuous and discrete NLS equations 

In this subsection we concentrate on the scalar NLS case, i.e. m = 1. Without restriction of 
generality, we can assume that the n x n matrix S has Jordan normal form. 

Lemma 6.5. 

(1) Let S be a single k x k Jordan block. For any k-component column vector V, there is a k x k 
matrix A that commutes with S and satisfies 

V = AVq where Vj = (0, . . . , 0, 1) . (6.33) 

A can be chosen invertible if the last component of V is different from zero. 

Furthermore, for any k-component row vector U , there is a kxk matrix B that commutes with 
S and satisfies 

U = UoB* where [/q = (1, 0, . . . , 0) . (6.34) 
B can be chosen invertible if the first component of U is different from zero. 

(2) Let S = block-diag{Si, . . . , Sr) with Jordan blocks Si. For any correspondingly decomposed 
n-component vector V" = (Vi , . . . ,Vj), there is an n x n matrix A that commutes with S and 
satisfies 

V = AVo with V,^ = {Vo]„ . . . ,V,]^) , (6.35) 
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where each vector Vq^j (of same size as Vi) has a 1 in the last entry and zeros otherwise. A can be 
chosen invertible if the last component of each Vi is different from zero. 

For any correspondingly decomposed n-component row vector U = {Ui, . . . , Ur) there is an nxn 
matrix B that commutes with S and satisfies 



U = UoB* 



with 



Uo = {Uo,i, . . • , Uo,r) , 



(6.36) 



where each row vector C/o,i (of same size as Ui) has al in the first entry. B can he chosen invertible 
if the first component of each Ui is different from zero. 

Proof: (1) A commutes with the Jordan block 5 if it is an upper-triangular Toeplitz matrix. Hence 



/ s 




s 



\ 







•• 

■•. 1 

s 



A 



( a\ a2 as 

■•• ■•. 



V 



\ 



■• as 

'■• a2 

ai / 



with arbitrary constants ai, . . . ,afc. It is now immediately verified that, for any given vector V, 
there is a matrix A of the above form such that AVq = V holds. A is invertible iff ai ^ 0, and 
ai coincides with the last component of V. Let ui, . . . ,Uk be the components of U and B the 
Toeplitz matrix built with the complex conjugates of these constants. Then U = UqB*. Clearly, B 
is invertible iff ui ^ 0. (2) immediately follows from (1). □ 

As a consequence of the lemma, we can always achieve that 



U = UoB*, V = AVo, 



(6.37) 



with matrices A, B commuting with S. As shown in section 6.1.1, the condition in the above lemma 
achieving that A and B are invertible is no restriction of generality, provided that S satisfies 
the spectrum condition (6.13). The latter implies that S is invertible. By a reparametrisation 
transformation^^ (and redefinitions of K and V in the discrete cases) the above family of solutions 
can then be expressed as 



Uq{I - e'^*K*'EK)-^E* S*-Wq, p = eUoE* K* {I - eEKE* K*)'^ E S'IVq , (6.38) 



with 



NLS 


V 


Sylvester equation 




continuous 





SK + KS* = VqUo 




semi-discrete 


1 


K - SKS* = VqUo 


C 5'^+^ g-ii (S+5-1-2/) 


fully discrete 


1 


K - SKS* = VoUo 


c 5^+1 [(/ + i (/ - s~^)){i - i (/ - s))-^Y 



where I = In and the matrices A, B only appear through the combination 

C = BA. 



(6.39) 



^*Tlie complex conjugation reduction requires A = A* and B = B* \n (5.19). 
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With the restriction to the focusing case, i.e. e = — 1, we now consider the special case where 
S is diagonal, hence 

S = dmg{si,...,sn), Uo = Vq' = {1,...,1) . (6.40) 

In the following we assume the spectrum condition (6.13), i.e. Sj + 7^ in the continuous case, 
and SiSj 7^ 1 in the discrete cases, for all i,j = 1, . . . ,n. According to proposition 6.4, with no 
restriction of generality we can then assume that the eigenvalues Si are pairwise different. As a 
consequence, A and B are diagonal, hence C is diagonal. Writing 

H = diag(Ci,...,^n), C = diag(ci, . . . ,c„) , K = (kij) , (6.41) 

we obtain 



NLS 


6 


kij 


continuous 






semi-discrete (AL) 




1/(1 - Sr.S*) 


fully discrete 


cs^' [(i+i(i-sr'))(i-i(i-so 




1/(1 - SiS*) 



K is a Cauchy-like matrix and Hermitian. ° 

Lemma 6.6. Let S be diagonal with pairwise different eigenvalues Si, subject to (6.13). 

(i) The matrix K of the continuous NLS case is positive definite if Re(sj) > for i = 1, . . . ,n 

(which means that S is positive stable). 

(a) The matrix K of the discrete NLS cases is positive definite if all Si lie in the open unit disk 
(which means that S is stable with respect to the unit circle). 

Proof: In case (i), the matrix K \s & positive semidefinite Cauchy matrix, see e.g. [146,149]. If the 
Si are pairwise different, Cauchy's determinant formula (see e.g. [146]) shows that K is invertible. 
Hence K is positive definite. In case (ii), we can write K = DCD^ with the Cauchy matrix 
C = (l/(Ai + A*)) and D = diag((Ai + l)/^/2, . . . , (A„ + 1)/V2), where Si = (A^ - l)/(Ai + 1) defines 
a bijection between the open unit disk and the open right half plane. Since the Si are assumed to 
be pairwise different, it follows that also the Aj are pairwise different. Hence C is positive definite, 
which then also holds for K. □ 

Remark 6.7. Since a diagonal matrix S is symmetric, the Sylvester equation in the continuous 
NLS case is the Lyapunov equation (6.16). The result formulated in part (i) of the preceding lemma 
can thus be deduced from a more general result in [148], which we recalled in remark 6.3. □ 

Proposition 6.8. Let S = diag(,si, . . . , s,i) with pairwise different Si, and Re(sj) > m the 
continuous NLS case, \si\ < 1 in the discrete NLS cases. Then the associated solution (6.38) of 
the (continuous, semi- discrete, respectively fully discrete) focusing NLS equations is regular (for all 
x,t £ R). 

Proof: Using a diagonal square root H^/^ of the diagonal matrix H, we introduce 

L = Hi/2A'(Hi/2)t. 

^^We actually arranged the form of the Sylvester equations in order to achieve this (at the price of introducing 77). 
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Since K is Hermitian and positive definite by the preceding lemma, also L is Hermitian, i.e. L'^ = L, 
and positive definite. Hence L possesses a positive definite Hermitian square root LV2 = (ii/2)t^ 
and we can write 

= {E'/^L'/^y [(^i/2)*^i/2]t^ 

For the expression in (6.38) that has to be inverted, this implies 

/ + {EKYiEK) = lV2)* (^/ + [(^1/2)^1/2] [(LV2)*LV2]t) (^1/2 Li/2y-i . (6.42) 

Since any matrix of the form / + AA^^ is invertible, this proves our assertion. □ 

Because of the following reason we can relax the assumptions in proposition 6.8 to the require- 
ments that the eigenvalues Si of the diagonal matrix S are pairwise different, and Sj + s* 7^ in 
the continuous case, SiS* 7^ 1 in the discrete cases. If we have matrix data {S,C,Uo,Vo), where 
S has an eigenvalue with negative real part (respectively which lies outside the unit circle), then 
we can apply a reflection symmetry and obtain new matrix data {S' ,U' ,V') that determine the 
same solution, and S' is given by S where the respective eigenvalue now has positive real part, 
respectively lies inside the unit circle. Furthermore, we can choose the data such that U' = UqB'* 
and V = A'Vq with invertible matrices A',B' that commute with S' . A reparametrisation trans- 
formation then achieves the form (6.38) with matrix data (5", C , Uq, Vq). Hence the assumption in 
the proposition that the eigenvalues of S have positive real parts, respectively lie in the open unit 
disk, is no restriction. 

6.3 Solutions of the matrix NLS equations 

By redefinitions of K and V in the discrete NLS cases, the solutions determined by proposition 6.1 
can be expressed as^^ 

q = U Z^^ S*~''V* , p = eUE*K* Z*~^ S~iV , (6.43) 

where 

Z = E*"'^ - e K*EK (6.44) 

and 



NLS 


7] 


Sylvester equation 




continuous 





SK + KS* = VU 




semi-discrete 


1 


K - SKS* = VU 




fully discrete 


1 


K - SKS* = VU 


[(/ + i (/ - S'^w - i (/ - s)r^Y 



This and what follows includes the case m = 1, i.e. solutions of the scalar NLS equations. For 
m > 1, the results are new. 



'If q and p solves the respective NLS system, then also —q and p. 
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6.3.1 Rank one single solitons 

Choosing n = 1, [/ is a column and V a row vector (with m components). Writing 

S = s, E = , U = u, V = v'^ , K = Kv'^u, 
with (m-component) column vectors u,v, we obtain^^ 



uv 



UV I 



where pq = — s~'^ uv'^ / (v'^ u) and 

NLS 7] 



P 



+ P0 



(6.45) 



(6.46) 



continuous 



semi-discrete 



fully discrete 







1 



1 



s+s* 



1—ss* 



1—ss* 



KV'Ue 



-s x—i s^t 



KV'US 



x+1 ~i{s+s-^~2)t 



+1 f 

l-i(l-s) 



Of course, we have to assume Re(,s) ^ 0, respectively \s\ ^ 1. In the focusing case, i.e. e = — 1, 
these solutions are regular and describe a single soliton with an attached "polarisation" determined 
by the vectors u and v. The solutions can be expressed as follows. 



q = csech{a {x — pt) — 6) e 
2 c uv^ 



i (px—ut- 



UV 



p 



+ P0 , 



where 



a + i/3 , c = a . 



2/3 , u = (3^ 



a 



in the continuous NLS case, with a, (3,6,6 G M. In the discrete NLS cases, 
s = e-"-''^, c = sinh(a)e^'5, -yT^, = ^-i g^+^+i (e-/?) 



and 



a 



sinh(a) sin(/3) , = 2 (1 — cosh(a) cos(/3)) 



(6.47) 



(6.48) 



(6.49) 



(6.50) 



in the semi-discrete case, whereas in the fully discrete case p and u are given in terms of a, /3 via 



l + i(l 



l-i(l 



-«-i/3^ 



(6.51) 



Since here q and p are rank one matrices, these are the most elementary matrix soliton solutions. 
Such a soliton has velocity p. Whereas a scalar (i.e. m = 1) NLS soliton carries a phase factor 
that can change via scattering with other solitons, here the phase is extended to the polarisation 
matrix. We shall see below how the latter changes in a scattering process. 



''Note that V = uv'' /{v''u) is a projector. Furthermore, V = uv'' / {v'' u)* satisfies VV — V. 
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6.3.2 Asymptotics of superpositions 

The following is a preparation for the elaboration of the asymptotics of multi-soliton solutions 
in section 6.3.3. We note that (6.25) can be written as Z = A2'^*~^B2 + AI'BBi with H = 
block-diagonal(Hi, Hg""*^) and 

(where / stands for the identity matrix of the respective size). Assuming that Ki and K2 are 
invertible, this implies 



B^^E-Ul -eK*-^E*-^K-^E-^) ^ A*-^ (6.53) 



where 



(cf. (6.32)), with the Schur complement of K2 (with respect to K), 

Li = Ki- Ki2K^^K2i . (6.55) 
In a limit in which — t- and H2 ~ H2 (asymptotic form), we obtain 



with the Schur complement of Ki (with respect to K), 

L2 = K2- K2iK^^Ki2 . (6.57) 

This allows us to compute the asymptotic form of q, given by (6.43), in this limit. In order to 
compute the asymptotic form of p, we write 

E*K*Z*-^ = ^ {BiE*K*Bl-^E*~^) (^I - e K'^ E'^ K*~^ E*-^^ ^ ylj-^ . (6.58) 

Since 

B,E*K*Br'E*-^ = ( +„5;,f p\5:r*"' """if^"^ ) , (6.59) 

we find 



In the same way, in a limit in which ^ — )• and Hi ~ Hi, we obtain 

Z-^ ~ B^i(^m{l-^LlE,L,Eiy' 0^^._i^ 

E*K*Z*-' ^ ^-1 / e^r' (/-e^iHJLJHi) ^ \ A^' . (6.61) 

V / 
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6.3.3 Superposition of rank one solitons 

Now we consider the focusing case, i.e. e = —1. We choose S = diag(si, . . . , s„) and write 



U={ui ■■■ Ur,) , V=\ : I , (6.62) 

with m-component column vectors Ui,Vi, i = 1, . . . ,n, and H = diag(^i, . . . , ^n). The corresponding 
Sylvester equation is solved by 

{ J'^i continuous NLS 
iJg.l* (semi- and fully) discrete NLS 

The NLS solution can now be expressed as 

n n 

q=^^niZE*r%s*-'^u,v], p= ^ {K~'),k[{Z*Ey%js-'^u,v]+po, (6.64) 

i,j=l i,k,j=l 

with the constant term pQ = —UK~^S~^V, dropped in the following, and 

n 

{ZE*)ij = Si, + k*, hj ^kCj ■ (6.65) 
fc=i 

For each Sj we introduce real constants Oj, f3j, pj,i>j as in (6.48), (6.49), (6.50), (6.51), and 

{si + s*)~^ continuous NLS , . 

(l-s,s*)-i discrete NLS ' ^ ' 

Furthermore, we choose^^ ai > 0, i = l,...,n, and^^ pi < p2 < ■ ■ ■ < Pn- Asymptotically, as 
t —7- — oo, the corresponding incoming solitons thus appear in decreasing numerical order along 
the X-axis. The outgoing solitons (as t — )• oo) then appear in increasing numerical order. To 
extract these asymptotic solitons from the solution in the limits t — )• ±00, we introduce a comoving 
coordinate Xi = x — pit, with respect to which the i-th. soliton is stationary up to a phase factor. We 
consider the limits t — t- ±00 while keeping Xi constant, for some fixed i. Then = e~"i(^~^'j*)+^j = 
g-a,[x,-(p,-p,)i]+5,^ follows that 



lim ^, = if J > i , lim = if j < i , (6.67) 

and 

lim = \i j >i, lim ^, = \i j < i . (6.68) 



In order to find the corresponding asymptotic forms of q and p, we express the solution as a 
superposition (6.22) with Si = diag(si, . . . , Si-i) and ^2 = diag(si, . . . , s„). Let Ui, U2 and Vi, V2 



^^This is not really a restriction since for a; 7^ it can be achieved generically by a reflection symmetry. 
^^For simplicity, we exclude the case where some of the pi are equal. 
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be the corresponding parts of U and V, respectively. Let Ki, K2 and Ki2,K2i be the corresponding 
diagonal and off-diagonal blocks of K. Using (6.56), (6.60), (6.61), with (^2)^ = Sik^u, we obtain 

' l + l7"f.P ^ (7^-^)*' l + ht'i^^' 7l-W "^x. = const. ^^'^^^ 



where we introduced 



— — '^i 



(6.70) 



and 



7 



(-) 



The latter vectors can also be expressed as follows. 



/ k 



i-l 



j,k=l 



det{Ki 



det 



11 



,(-) 



i,fc=i 



^j— 1,1 ■ ■ ■ ^i— 1, 
\ Itl • • • 

/ ^11 

det 



.(6.72) 



\ kii ■■■ fcj,j_i s- 'Vi j 

The first equation e.g. can be verified by Laplace expansion of the determinant of the big matrix on 
the right hand side with respect to the last row, and subsequent Laplace expansion of the cofactors 
with respect to the last column. The expressions in (6.72) are examples of quasideterminants (see 
e.g. [150]). Note also that 

/fell ••• kxi \ 



{L2)ii = {K2-K2iK^^Ki2)i 



1 



det 



(6.73) 



det(/^i) 

\ ki\ ■ ■ ■ kii / 

so that all constituents of the formulae in (6.69) have nice expressions in terms of the vectors Ui 
and Vi, and kji with j,l = 1, . . . ,i. The expressions in (6.69) are arranged in the form of the single 
soliton solution (6.46). This is more clearly seen by noting the following result. 

Proposition 6.9. 

v(-^^u(-' = 7/-^ vj Ui . (6.74) 
Proof: In the continuous NLS case, we have vjuj = {si + s*) kij. Hence 

i-l i-l 

j,k=l l,r=l 
i—1 i—l i—1 

= vjui- ^vJui{K{'^)irkri- vlUikij{K{^)jk+ ^ vluikij{K{'^)jk{K{'^)irkri 

l,r=l j,k=l j,k,l,r=l 

i—l i—1 i—1 

= (si + s*) (hi - ^ kij{K^^)jkkki] - ^ kij{K{'^)jkSk [kki - ^ kki{K{^)irkr 



j,k=l jM=l 
i-l i-l 

^ [kii - X] f^iji^l^)jkhl^Sj: {Ki'^)lrkr 
l,r=l j,k=l 



Lr=l 
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Now we observe that the last two sums vanish as a consequence of the identities 



kri ^ ^ krj{K^ )jkkki — 0, kir ^ ^ kij{K^ )jkkkr — ^OT r <. i . 
j,k=l j,k=l 

The first identity is e.g. obtained from the first identity in (6.72) by replacing Uj by krj, j 
Hence 

(-)T (-) {^2)ii T (-) T 

In the discrete NLS case, using vjuj = (1 — sis*) kij, we obtain 



i-l 



i-l 



(1 - SiS*) {L2)ii - Si ^ kij{K^ ^)jk ^ [kki - ^ kki{K^ ^)irkr 

j,k=l l,r=l 

i-l i-l 

-Si ^ (^kii- ^ kij{K^^)jkkki^Si {K{^)irkr 

l,r=l j,k=l 



Now we can rewrite (6.69) as follows, 



In the same way we obtain 



P ■ 



»,{-) ^(-)T 
-V i i 



as 



t —7- — oo 
Xi = const. 



(+) „,{+)t 



p ■ 



1 + hrc 



(+)A|2 -y(+)T^(+) 



as 



t +00 
= const. 



where 



(+) 



(^i) 



II ( + ) 



kit 



The vectors can also be expressed as follows. 



(+) 



u 



det(A'2 



det(A:2) 



det 



det 



/ Ui 

ki+\^i 



Un \ 

ki+l.n 



\ kni ■ ■ ■ kfiYi J 
^ ^i ^"^i ^i,i+l ■ ■ ■ kin ^ 

\ Sn^Vn A;n,i+1 " " " ^nn / 



and satisfy 
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which is proved in the same way as proposition 6.9. Furthermore, we have 



det{K2) 



det 



(6.80) 



knn j 



All this nicely expresses the constituents of the formulae (6.76) in terms of the vectors Ui and i^j, 
and kji with j,l = i, . . . ,n. 

Clearly, if i = 1 respectively i = n, we have 



1^1 it *—r] '^1 ■'^1 

Si 



K 



-1 



1 + 161 
^-1 



-r, Ui V{ 
'1 

v{ui 



as 



1 + 



- 2 



t — )• — OO 

xi = const. 

as _ , (6.81) 

Xr, = const. 



By comparison of the asymptotic expressions (6.75) and (6.76), the solitons experience a position 
shift and a change of polarisation through the scattering process. 

6.3.4 2-soliton solutions 

If 71 = 2, we have 



S 



si 

S2 



u 



Ui U2 ) 



V 



(6.82) 



Using a reparametrisation transformation (with diagonal transformation matrices) and a transfor- 
mation (5.20), we can achieve that 



{vjuj)* = v]Ui 



(6.83) 



so that K is Hermitian. But now we have to admit in the expression for a factor Cj G C. The 
matrix K, given by (6.63) in terms of the eigenvalues Sj and the vectors Uj,Vj, j = 1,2, has the 
form 



K 



ki k 
k* k2 



(6.84) 



(where the ki are real), and we obtain the following explicit solution. 



1 r 



^ = D 



;i + H i^r + ^'6^2) eist'u.vi + (1 + kf leii^ + f^^*^) ^2*4"V^^ 



-{k fci ei + k* k2 6) 6*^2* (st^U^vl + Sl-'^U2VI] 



D 



ki + det(K) k2 |6r) ki\ii\^s^'^uiv\ + [k2 + det(K) ki A^al^r^a ""2^ 



+ {k*-k det(K) 6^2*) eii2S2''uivl + {k-k* det{K) Cl^2) CiC2S^''u2vl 



where 



D 



k* 



(6.85) 



(6.86) 
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Remark 6.10. If 5" = s/„ with n < m, we have a "rank n sohton". It follows from proposition 6.4 
that we can choose all m-component vectors composing U, respectively V, linearly independent. 
In particular, this means that there are at most different such solitons. If n = m, then with 
a reparametrisation transformation we can achieve that U = V = Im, while replacing H by CH 
with a constant m x m matrix C. Imposing the condition (6.83) on the vectors Ui,Vj, so that K 
is Hermitian, and choosing C symmetric, the resulting NLS solution is regular. In the continuous 
NLS case, q is then symmetric. Superposing such rank m solitons (with symmetric Cj), we obtain 
again continuous m x m matrix NLS solutions with symmetric q. This may be of interest for the 
aforementioned spin 1 Bose-Einstein condensate model. □ 



6.3.5 Regularity condition for matrix solitons in the focusing NLS case 

We address the question whether a superposition (in the above sense of "superposing" the cor- 
responding matrix data) of regular solutions again leads to a regular solution. Hence we should 
check the determinant of (6.25) for the occurence of zeros. As an intermediate step, we compute 
the Schur complement of the first block-diagonal entry of Z, which is 

2) := ^2 + ^2*1^1^12 - + Kl{EiKi){Zi + Kl2^2K2i)-\KlEiKi2 + K^2^2K2) .(6.87) 

Here we assume that, for the superposed data, a matrix K as in (6.22) exists such that (6.23), 
respectively (6.24) holds. A sufficient condition is (6.13). We further assume that Zi and Z2 are 
(for all X, t) invertible, so that the constituents of the superposition are regular, and also that 
+K2iZ^ is invertible. By the matrix inversion lemma (see e.g. [137], Corollary 2.8.8), the 
latter condition guarantees the existence of 

(Zi + Ki*2S2i^2i)"' = - 1^1*2(^2 ' + K2iZ^^Kl2r^K2iZ^^ . (6.88) 
Using this formula in the preceding equation, and noting that 

(^2 ^ + K2lZ^'^Kl2)~^ K2lZ{^Ki2'E2 = "2 — ("2 ^ + K2lZ{^Ki2)~^ 

= E2K2iZ{'K*,2{E^' + K2iZ^'Kl2r\ (6-89) 

we obtain 

D = Z2 + K^iEiKi2-K;^EiKiZ^^KlEiKi2 

+K;,EiKiZ{'KUE^' + K2iZ^'K^2r'K2iZ^'KlEiKi2 
-K*2{E^^ + K*2iZ{'K*,2r'K2iZ^'KlEiKi2 
-K2*iSiKiZf 1^1*2(^2 ' + K2iZ^^Kl2r^K2 

-K;E2K2 + K^{E^^ + K2iZ^^Kl2r^K2 . (6.90) 
Using the identity EiKiZ^^K^ = I — Z^^^E^^ in the third summand, this can rewritten as follows, 

2) = {K^^EiKiZ^^Kl2 - K;KE^^ + K2iZ^^Kl2)-HK2iZ^^KlEiKi2 - K2) 

+ Kl^Zl~^Ki2 . (6.91) 

With the help of the identity EiKiZ^^ = {Z^^KlEi)* ^ we obtain 

S = {K2iZ^^KlEiKi2 - i^2)*(H2 ' + K2iZ^^Kl2r\K2iZ^^KlEiKi2 - K2) 

+(S2-^ + K2iZ^^Kl2r . (6.92) 
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If the data {S2, U2, V2) belong to a 1-soliton solution, writing ^2 = 5,^2 = ^, U2 = u, V2 = v'^ , 
K2 = k, K12 = ki, and K21 = ^2' have 

f Zi+ik\kl KlEiki+ikk\\ 

\^ k*kl + fc^HiA'! Z2 + fc^Hifci ) ' ^^-^""^ 
where Z2 = + Using determinant identities (see e.g. [137], Fact 2.14.2), we find 

det(Z) = det{Zi + ^klkl)^ = det{Zi){l + ^klZ^^kl)^ 

= C det(Zi) (r' + klZ^'kl) ((r' + kiz^'kir 

+{klZ^^KlEiki - k)* {C^ + klZ^^k\)-^{klZ^^Kl'~iki - A:)^ 
= i det(Zi) + klZ{^k\\^ + \klZ{^KlEiki - k\A . (6.94) 



This can only vanish if /c = k\Z^^ K^r^iki and simultaneously = —k2Z^^k\ at some value of 
X and t. Generically it is not possible to satisfy both equations at any value of x, t. (Note that 
these equations consist of four real equations.) In particular, this shows that the matrix n-soliton 
solutions are generically regular. 

6.3.6 Cases in which the Sylvester equation has not a unique solution 

Several of our results so far in this section rest upon the assumption that the condition (6.13) is 
satisfied, in which case the Sylvester equation in proposition 6.1 possesses a unique solution. Now 
we address those cases in which this condition is violated, confining ourselves to the continuous 
matrix NLS equation (6.3). The simplest such cases occur if the matrix S consists of a single 
Jordan block with a purely imaginary eigenvalue. By application of a Galilean transformation (see 
remark 6.2), this case can be reduced to that where the eigenvalue is zero. 

Example 6.11. (1) Let n = 1 and 5 = 0. Writing U = u and V = v'^ with m-component vectors 
u, V, the Sylvester equation requires v^u = and leaves K = k completely arbitrary. If m = 1 (the 
scalar NLS case), the corresponding NLS solution vanishes identically. But if m > 1 (the matrix 
NLS case), this is not so^*^ and we obtain the constant rank one matrix NLS solution 

^——:ruv'^ where v^u = Q. (6.95) 



A Galilean transformation generalises it to 



1 ei/3{^-/3t)^^t^ (6.96) 



' l-e|A;| 

with /3 G M. This is what we would have obtained had we started above with 5" = i /3, see remark 6.2. 
As a consequence of v^u = 0, we have qq* = 0, so that q solves the linear and the nonlinear parts 
of the NLS equation separately. 

(2) Let n = 2 and S a Jordan block with zero eigenvalue, hence 

5=(0 M, U={m a2), V=(ll)- (6.97) 



^''This is in contrast to the Hermitian reduction case treated in section 7, where the corresponding Sylvester 
equation also leads to the trivial solution if m > 1. See remark 7.2. 
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Figure 1: The norm \\q\\ of the 2x2 matrix focusing NLS solution, determined by the data (6.99), 
at t = —4,-2,0,2,4. This is a rational "matrix breather". 



The Sylvester equation only has a solution if &2 =0 and &| ai = &2 ^2- Let us more concretely 
consider the case m = 2. The solution K of the Sylvester equation depends on two arbitrary complex 
constants (expressing the non-uniqueness of K). The resulting components of q are quotients of 
a linear polynomial in x by a quadratic polynomial, hence q vanishes as — t- oo. Applying a 
Galilean transformation (see remark 6.2), these solutions become t-dependent.^^ A special regular 
member from this family is 



q 



l + {x-2l3t) 



1 

- 2pt 



-{x-2/3t) 
1 



/3e 



(6.98) 



(3) For n > 2, with 5 still consisting of a Jordan block with eigenvalue zero, one obtains t-dependent 
rational solutions. This includes regular solutions with a breather-like behaviour. For n = 4, a 
special solution of the 2x2 matrix focusing NLS equation is determined by 



S 



/ 1 \ 

10 

1 

\ / 



u 



1 

1 



V 



1 ^ 






( 
















-1 


, K = 




-1 



















-1 


1 








V 1 


^ ) 




\ 





1 


-1 


/ 



This rational solution is regular, vanishes as 
Fig. 1. 



X 



oo. 



(6.99) 

and shows a breather-like behaviour, see 

□ 



Of course, we can superpose any number of solutions of the kind described in the last example, 
see the next example, and also compose such data with (e.g. solitonic) data satisfying (6.13), as 
in the next but one example. In doing so, we have to take the parameters /3j into account, since 
a Galilean transformation applied to the composite solution can only eliminate or create one of 
these parameters. In the special case where S = diag(i /3i, . . . , i/3„), with real (3i ^ 0, the resulting 
solutions are periodic in x and t, since they are rational expressions in trigonometric functions. 

Example 6.12. The following superposition of data of the kind considered in the preceding ex- 
ample determines a regular solution that vanishes as |x| — )• oo, see Fig. 2. 



S 



/ if3 1 \ 

i/3 

1 

\ / 



u 



10 10 
10 1 



, V 



1 1 


o\ 









i 
2 







1 




1 








i 


1 


-1 


i 
2 


2 





2 




Vo 


1 / 




^ 


i 
2 


1 


-1 / 



□ 



^The absolute values of the components are still rational functions, hence these solutions resemble the "rational" 
soliton solutions known in the scalar NLS case [151]. 
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Figure 2: The real part of qu for the matrix NLS solution q determined by the data in example 6.12 
at t = -4, ... ,4. 

Example 6.13. Let S = diag(s,0), s iM, and U,V as in (6.97) with m = 2. The Sylvester 
equation then admits a solution K iff &2 ^2 = 0. Writing aj = (aji,Oj2), bj = {bii,bi2), we solve 
this constraint by setting 621 = ca22 and 622 = —ca2i with c G C. Then K takes the form 



K 



aiifoil+ai2fei2 

s+s* 
aiia22— ai2a2i 



a2ifcii+a22fei2 
s 

k 



(6.100) 



where G C is another arbitrary constant. The resulting solution family of the 2x2 matrix 
focusing NLS equation depends on nine complex parameters and one real. Fig. 3 shows a regular 
solution from this family. □ 

Recall that any NLS solution corresponding to matrix data (5, U, V) consisting of several blocks 
can be regarded as a (nonlinear) superposition of solutions corresponding to these blocks, say 
(S'i,C/i,Vi) and {S2,U2,V2)- If '^('S'l) H (t(— S*!) = in the continuous NLS case (respectively 
a{Si)na{S2^) = in the discrete cases), the solution is completely determined by the "constituent 
solutions". Soliton solutions are the examples par excellence. If cr(5i) n o"(— 5^) / (respectively 
a{Sl)na{S^^) / 0), there is no unique solution for the off-diagonal blocks of the matrix K solving 
the Sylvester equation for the superposition. The parts then no longer determine the composite 
solution completely, since additional freedom enters through arbitrary parameters in K. But if 
the eigenvalue of Si has a non-zero real part in the continuous NLS case, respectively a modulus 
different from 1 in the discrete NLS cases, then we can apply a reflection symmetry to S2 (provided 
that the associated matrix K2 is invertible). This leads to equivalent data involving Jordan blocks 
Si and 5*2, where now (t{Si) n (t(— Sg*) = 0. This means that such a NLS solution can be realised 
in different ways as a superposition, i.e. with different constituents. 

7 Hermitian conjugation reduction of the continuous NLS system 

By decomposition, the Hermitian conjugation reduction condition (3.33) for the continuous matrix 
NLS system becomes 



eq' 



p = p' 



P = P 



(7.1) 
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Figure 3: A sequence of snapshots of a solution from the family in example 6.13. The plots 
show the absolute values of the components of the 2x2 matrix focusing NLS solution q at times 
t = —.2, 0, .2, .4, .6. Here we chose the parameter values s = 2 (1 + i), c = A; = an = 022 = 1 and 
021 = -1, ^12 = i, ai2 = bii = 0. 

where e = ±1. It reduces the NLS system to the mi x m2 matrix NLS equation (see e.g. [3]) 

igt + - 2egg'l'g = 0, (7.2) 

supplemented by px = —eqq^. 

If 1712 = 1 and mi > 1, (7.2) is the vector NLS equation (or coupled NLS equations) 

iqt + qxx-2e\q\^q = 0, (7.3) 

first treated in [77]. 

If mi = 7712 and if q is symmetric or skew-symmetric, i.e. q^ = ±q, then (7.2) coincides with 
(6.3). This includes the special case of a spin-1 Bose-Einstein condensate mentioned earlier. For 
the skew-symmetric case, see e.g. [152]. 

Proposition 7.1. Let S, V, V be constant matrices of size n x n, nx mi and n x m2, respectively, 
and K, K constant Hermitian solutions of 

SK + KS^ = VV\ S^K + KS = VVK (7.4) 

Then 

q = V^ (e^-^ -eKEKy^ V where E = e"^ ^''^ "^^ (7.5) 

solves the mi x m2 matrix NLS equation (7.2). □ 

Proof: Setting ni = n2, renamed to n, we recall from section 4.4 the reduction conditions to be 
imposed on the matrices entering the solution formulae. With the block structure expressed in 
(4.23), (4.48) and (5.9), the statement then follows directly from proposition 5.1. □ 
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The class of matrix NLS solutions provided by proposition 7.1 has also been obtained in [98, 153] 
via inverse scattering, including conditions under which the solutions are everywhere regular and 
exponentially localised (see also [143]). Multiple soliton solutions of the matrix (including vector) 
NLS equation appeared previously also in [3,77,81,82,87-94,117], for example. 

If K and K solve the Lyapunov equations (7.4), then also K'^ and K'^ . As a consequence, there 
are then also Hermitian solutions of these equations. We further note that the equations (7.4) have 
unique solutions if 

a(5)na(-St) = 0, (7.6) 
in which case K and K are then necessarily Hermitian. 

Remark 7.2. This concerns cases in which (7.6) does not hold. If 5 = 0, then (7.4) implies 
F = V" = 0, leading to g = 0. If 5 is an r x r Jordan block with eigenvalue zero, i.e. S., 



more involved inspection of (7.4) shows that K is left-upper-triangular, K right-lower-triangular, 
= (I'l, . . . , Vk, 0, . . . , 0) and = (0, . . . , 0, t);, . . . , 1?^), where k = r/2, / = A; -|- 1 if r is even, 
and /c = (r — l)/2, / = /c -|- 2 if r is odd. (7.5) again leads to g = 0. Corresponding solutions of the 
complex conjugation invariant matrix NLS equation, obtained in section 6.3.6, have therefore no 
analogue in the Hermitian invariant NLS case. □ 

Remark 7.3. The first result noted in remark 6.3 for symmetric S can also be obtained from 
proposition 7.1 by setting V = ^V* and K = . □ 

Example 7.4. Let mi = m2 = 1, hence we consider solutions of the scalar NLS equation. Let 
S = diag(si, . . . , Sn) with pairwise different eigenvalues Sj. According to lemma 6.5, without 
restriction of generality we can set 



V = A 



( 1 



V 1 



V = A 



( 1 \ 



V 1 / 



(7.7) 



with matrices A^A satisfying [^,5"] = = [^, /S^]. Using an argument similar to that in sec- 
tion 6.1.1, A and A can be chosen invertible. By application of a reparametrisation transformation^^ 
the solutions can be expressed as 

. / 1 



1 



)(= 



■t-i 



\ 1 



where 



and K, K have to solve 



E = Ce 



with C = A^A , 



/ 1 



SK + KS^ 



1 \ 



K + KS 



(7.; 



(7.9) 



(7.10) 



\l ... 1/ 

Since 5 is symmetric (i.e. S'^ = S) and K = K'^ Hermitian, the solutions coincide with those ob- 
tained via the complex conjugation reduction in section 6, yielding the well-known soliton solutions 
of the scalar focusing NLS equation. □ 



The Hermitian conjugation reduction requires B = At and B = A'^ in (5.19). 
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Remark 7.5. The Hermitian conjugation reduction does not work for the discrete matrix NLS 
systems, except when mi = m2 and qq^ = q^q = |xP-^ with a scalar x [3,92, 118, 123, 124, 130] 
(see also [119,120]). One is then led to the semi-discrete vector NLS equation (or coupled NLS 
equations) 

i^+q+ -2q + q- -e\q\'^{q+ +q-) =0, (7.11) 

where the vector q is obtained via expansion of the matrix q in terms of a basis of an appropriate 
Clifford algebra. A corresponding analysis in our framework will be postponed to a separate work. 
In appendix D we present a non-local modification of the Hermitian conjugation reduction which 
does work for the semi-discrete NLS system. □ 



7.1 Solitons of the focusing matrix NLS equation 
7.1.1 Rank one single solitons 

Setting n = 1, proposition 7.1 determines rank one single soliton solutions of the focusing mi x m2 
matrix NLS equation. Writing V = and V = v\ with a reparametrisation we can achieve that 
the TTT-i-component column vector v and the m2-component column vector v satisfy llt^ll^ := v^v = 1 
and jji'll^ = 1, at the price of having to replace H in (7.5) by 

^ = ce-'^'-'''' = 2ai where ^ = e~'i^ -^')') e-''^^-^^'^+^ , (7.12) 

with c > and s = a + i/3 with a, /5 G M, a > 0.^^ Furthermore, we wrote c = 2ae^ with 5 € M. 
The solutions of the two equations (7.4) are then given by 

K = K = K where k = = — , (7-13) 

s + s* 2a 

and the matrix NLS solution (7.5), in the focusing case e = — 1, takes the form 

q = ^^^vv^ = asech{a{x-2pt) + 6) e' (/3^+("'-/3')<) . (7.14) 

This describes a soliton with amplitude given by a, with velocity 2/3 and polarisation matrix vv"^ . 



7.1.2 Superposing solutions 

Given two data sets {Si,Vi,Vi), i = 1,2, that determine solutions according to proposition 7.1, 
with corresponding Hermitian solutions Ki,Ki of the corresponding Sylvester equations, we can 
superpose them, 

where Ki2,K2i and ^12,-^^21 have to solve 

S1K12 + K12SI = ViV^ , S2K21 + K2iS\ = V2VI , 

S\Ki2 + K12S2 = ViV^ , SIK21 + ^21^1 = V2VI . (7.16) 



^^a 7^ is required in order to have a solution of the Sylvester equations (7.4). a > can then always be achieved 
by a reflection symmetry. 
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This system reduces to only two equations by setting 

K2i=k\^, K2i=k\^. (7.17) 

We have 

K"K=( _Zi - e K12E2K21 -eKiEiKu- eKi2E2K2 \ , . 

' ^ ~ \ -eK2E2K2i-eK2iEiKi - e ^21^1/^12 J' 

where Zi = - cKiEiKi. 

The fohowing paves the way towards the asymptotic structure of multi-sohton solutions, worked 
out in section 7.1.3, in analogy to the analysis in the case of the complex conjugation reduction in 
section 6.3. First we note that Z = AE'^~^A + BEB with E = block-diag(Hi, ^2"^) and 

Assuming that Ki and Ki are invertible, in a limit in which E^^ — )■ and H2 ~ H2, we obtain 

with the Schur complements of Ki (with respect to K) and Ki (with respect to K), 

L2 = K2- K2iK-^Ki2 , L2 = K2- K2iK-^K^2 , (7.21) 
which are also Hermitian. 

7.1.3 Superposition of rank one solitons 

A solution of the mi x m2 matrix NLS equation composed of n rank one solitons is obtained from 
proposition 7.1 by choosing 5 = diag(si, . . . , s„) with Sj = + i/3j. We assume > (which can 
be achieved generically by a reflection symmetry, see also [143]). Let us write 



(7.22) 



where Vi and Vi are mi- respectively m2-component column vectors. Conditions (7.4) are then 
solved by 

K = (%) , K = ikj) where % = , hj = • (7.23) 

Si + Sj Sj + Sj 

K and K are Hermitian. Using a reparametrisation transformation, we can achieve that the vectors 
Vi,Vi are normalised, i.e. 

Il^^ill = 1 = ll-^ill ■ (7.24) 
Hence ka = ku = 1/(2 Oj). Now E in (7.5) has to be replaced by H = diag(^i, . . . ,^n) with 
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where we set c, = 2aj e^^ . 

Next we choose Pi < P2 < ■ ■ ■ < Pn-"^^ As t — )• —00 the sohtons are then arranged along 
the X-axis in decreasing numerical order (recall that the velocity of the i-th. soliton, when isolated 
from its partners, is 2 Pi). After the interaction, as t — t- +00, the order is reversed. In terms of 
the i-th comoving coordinate Xi = x — 2 Pit, we have = 2a j exp(— Qj Xi + 2{Pj — Pi) t] + 5j). 
Hence, as t — >• —00 (keeping Xi constant), we have — )■ for j > i and — )• for j < i. 
Furthermore, as t — )• 00, we have — )• for j > i and £,j — )■ for j < i. In order to find 
the corresponding limit of q, for fixed i, we express the solution as a superposition (7.15) with 

= diag(si, . . . , Si-i) and S2 = diag(si, . . . , s^). Let Vi, V2 and Vi, V2 be the corresponding parts 
of V and V, respectively. Let Ki, K2 and Ki2,K2i be the corresponding diagonal and off-diagonal 
blocks of K, and correspondingly for K. Then we can use the results of section 7.1.2, in particular 
(7.20) with {E2)ki = Ci^ik^a, in order to work out the limits of q given by (7.5) as t — )• —00, 
respectively t — t- -|-oo. To formulate the results, we introduce new vectors via 



h, 7, 
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det{Ki 



det 
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( Ki ■ 
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det 



hi \ 
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(7.26) 



and 
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/ Vi 
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where the constants r], 
hold, 



^ \Vi 



T]-"''^ are chosen such that the following normalisation conditions 
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Proposition 7.6. 



n 



(-)\2 



ku' {L2)ii 
-- k^^" {L2)ii 



k^ {Li)ii 
-- k^^" {Li)ii 



(7.29) 



with the Schur complements Li 
(7.21).^^ 



Ki - K12K2 ^K2i, Li = Ki- K12K2 ^K2i, and ^2,^2 given by 



^*Here we leave aside the cases where some of the /3i coincide. 

^^Recall that Ki,Ki are Hermitian and (7.17) holds, hence Li, Li are Hermitian and the diagonal entries are 
therefore real. 
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Proof: Based on (7.23) and (7.24), the proof is similar to that of proposition 6.9. □ 
In particular, since we chose Oi > and thus ka > 0, this implies {La)ii > 0, {La)ii > 0, a = 1, 2. 

Let 



fi = \/{L2UL2h , ir = kr,' \/iL,UL,)u . (7.30) 
Then the following asymptotic relations hold,^^ 



1 + hrc 



(-)t.|2 



vl VI as t — )• — oo , Xi = const. , (7.31) 



q ^ a ^1 as t ^ +oo , Xi = const. . (7.32) 

These asymptotic expressions determine the changes in relative positions and polarisations of the 
solitons due to the collision. For the vector NLS equation, such results have already been obtained 
in [3,91,93] in a different way. 

7.1.4 2-soliton solutions 

For n = 2 we can write 

where 

'« = -^ = :^, ^ = 4^, J = 4^. (7.34) 

Si + s* 2ai si + si + 

The corresponding mi x m2 matrix focusing NLS solution is then explicitly given by 

^ \ei + «2i6i' Ci + kk i6i'C2) viv\ + (e; + ^i\ci\H*2 + k*k*\a\^ei) V2vi 



D L 



-{kKl\Clfe2 + k*K2\C2\^ei) Vivl - {kKi\ClfC2 + k*K2\C2\^ei) V2vl , (7.35) 

where 

'^l/. H-|2Mt |2\/'i , '^2. |F|2nia i2\ , 1 I , . 7_*a |2 



I)= l + ^(KiK2-|A:|")|6r 1 + — (^i^2-|A:r)|6r + '^iA:6 + '^2A:*6 • (7.36) 

As a consequence of the definitions (7.34) and oii > 0, we have kiK2 — |A;p > and KiK2 — \k\'^ > 0, so 
that the solution is regular. If the vectors i^i and V2 are orthogonal (with respect to the Hermitian 
inner product), and also vi and V2, then K and K are diagonal and the superposition is simply 
the sum q = qi + q2 oi the single soliton solutions. 



^Clearly, if i = 1 respectively i — n, keeping respectively constant, then we have 

q ~ — i-^-^ — Vi Vi as t — 7> — oo , q ~ = v„ v' as t — >■ +cxj . 
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For the above solution we have the fohowing asymptotic relations 



where = Ki ^i, 73 ' = 7[+' = ^/deUK)det{K) / {ki K2), and 
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8 Conclusions 

An appealing aspect of the bidifferential calculus approach is the simplicity of the basic structure 
one starts with, as compared with the typical complexity of a Lax pair. Our work demonstrates 
this in case of NLS systems and their discretisations (see also appendix E). The step from the 
continuum NLS bidifferential calculus to its discretisations is quite evident and the general solution 
generating method of section 2 is easily adapted. In this work we concentrated on the application 
of this method, which seems to be the most direct and simplest method to generate a large class 
of relevant exact solutions. In [1] we have shown how other methods, like Backlund and Dar- 
boux transformations, universally arise in the bidifferential calculus framework. The basic general 
result formulated in theorem 2.3 may actually be regarded as a shortcut to the binary Darboux 
transformation method [1]. 

In this work we presented infinite families of exact solutions of matrix continuous, semi- and 
fully discrete NLS systems and some of their reductions. The solutions are parametrised by con- 
stant matrices, and the only task left is to solve a Sylvester equation (which appeared as a "rank 
condition" in the literature). In the focusing case, they include the soliton solutions. But there 
are more and a complete understanding of the relations between matrix data and the behaviour of 
solutions has not yet been reached. But the expression of solutions in terms of matrix data turned 
out to be very convenient in this respect. The regularity proofs and the asymptotic analysis for 
some classes of solutions in sections 6 and 7 provide good examples. 

Besides a comparatively short treatment of the previously (with different methods) studied 
matrix continuous NLS equation with a Hermitian conjugation invariance, we considered in great 
detail another matrix NLS equation with a complex conjugation invariance. They coincide in 
the scalar case. But in the square matrix case, the latter has discrete counterparts, whereas 
the Hermitian conjugation invariant matrix NLS equation seems to admit a discretisation only if 
additional constraints are satisfied [3,118,123,124,130] (but see also appendix D for a corresponding 
non-local version). We therefore concentrated on an exploration of continuous and discrete matrix 
NLS solutions of the complex conjugation invariant NLS equations. We would also like to generate 
solutions of vector discrete NLS equations [3,118,123,124,130], using the framework presented here, 
and plan to elaborate on this in more detail in a separate work. Of further interest is the possibility 
to extend the present work in order to obtain solutions of (2 + l)-dimensional generalisations of 
matrix NLS equations. 

In the defocusing NLS case we have not yet reached, via the method of section 2, the important 
class of dark soliton solutions (see in particular [39,81,87,154-160]) and also rational solutions as 
e.g. obtained in [161]. A generalisation of this method seems to be required. 
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Appendix A: Solutions of the matrix KdV equation 

Let Bq be the algebra of smooth complex functions on E?, and B its extension by adjoining an 
operator dx that satisfies dxf = fx-\- f dx (where fx denotes the partial derivative of / with respect 
to x). We define linear maps d,d : B^B® via 

d/ = [dx, f] Ci + \[dl f] C2 , d/ = -\[dl f] Ci + l[dt- dl f] (2 , (A.l) 

where Ci)C2 is a basis of /\^(C^). In particular, we have d{dx) = d{dx) = 0. The maps d,d extend 
to linear maps B ® B A('C^) satisfying (2.3) and the graded Leibniz rule. Moreover, 

they extend to Mati(5) ® A(C^). For cp £ Mat(m, m, Bq), we have 



and thus 



del) = (l)x Ci + \{(t)xx + 2 cl)xdx) C2 , (A.2) 

d(/> A d(/> = -]^{(l)x^)x Ci A C2 , dd0 = {^\^(l^xt + -^(l^xxx^ Ci A C2 , (A.3) 

so that (2.14) turns out to be equivalent to the m x m matrix KdV equation 

1 3 

ut - ^Uxxx - -^{u'^)x = where u = (px ■ (A.4) 

In the following we apply results of section 2 in order to generate exact solutions of (A.4). 
Inspection of the linear equation dX = (dX) P suggests setting P = S — Idx, where d- and d- 
constancy requires that S does not depend on x and t. The linear equation is then equivalent 
to 

= — 2Xa:;S , Xt = Xa;S^ . (A. 5) 

These equations make sense for S G Mat(n, ?i,C) and X G Mat{n,n,BQ). They are solved by 

X = Ao + Ai H , H = e-2(-s+ts3) ^ (^.6) 

where Aj are constant complex matrices. Assuming Aq invertible, application of theorem 2.3 
reduces to an application of corollary 2.5. (2.18) takes the form 

SK + KS = VU, (A.7) 

and (2.20) becomes 

(/) = UH(I-KH)-^V. (A.8) 

This solves the mxm matrix KdV equation for any matrix data (S, U, V) for which (A.7) admits a 
solution K. Under the additional assumption that S is positive stable, so that e~^^ — t- as y — )■ oo, 
the condition^^ (A.7) is solved by [131,162] 



K 



/ e-y^Vte-y^ dy . (A.9) 
Jo 



^^For m = 1, (A.7) is a well-known rank one condition in the context of the scalar KdV equation. See [141, 142] 
for a corresponding general analysis in the framework of bounded linear operators on a Banach space. 
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In the scalar case, i.e. m = 1, we can apply more directly the formulae of theorem 2.3 to obtain 

(p = UYX-iV = tr(VUYX~i) =tr((XS + X^-RX)X"^) 

= tr(S-R)+tr(X^X-^), (A.IO) 

where R = H + dx- The second step requires that the matrices do not contain operator terms 
(involving dx)- Dropping the irrelevant constant summand, this becomes 

cP={logT)x, (A.ll) 

where 

r = det (l - e-^«Ke---«-2*s^) = det (l - e-^^VUe-^^e-^*^' dy^ . (A.12) 

In the last step we inserted (A. 9), assuming positive stable S. This determines the family of (here 
in general complex) KdV solutions derived in [134] (see (4.1) therein). 

Alternatively, using lemma 6.5, (A. 8) can be expressed in the scalar case (m = 1) as (j) = 
V5 C H (I - K C H)-i Vo, with V(5 = (1, . . . , 1). Here K now has to solve SK + KS = VoVj, and 
C is any n x n matrix over C that commutes with S. If S is diagonal and if the eigenvalues satisfy 
•Si + Sj ^ for all i,j = 1, . . . ,n, then K is the Cauchy matrix with entries l/{si + Sj). 

Appendix B: A relation between the semi-discrete NLS and modi- 
fied KdV equations 

Let us recall from (3.19) the m x m matrix equation for l/( in the semi-discrete NLS (AL) case, 

iJU -2U +U+ {I -U"^) + {I -U^)U~ = 0. (B.l) 

Substituting 

U = i{Uf e-^'-^^V xeZ, (B.2) 

with a new matrix variable V, we find28 

V + V+ {I + V^) - {I + V^)V- =0, (B.3) 

which is the m x m matrix (semi-) discrete mKdV equation. As a complex matrix equation, the 
latter is thus an equivalent form of (B.l). 

Proposition B.l. Let S,U, V be complex matrices of size n x n, m x n and n x m, respectively, 
and K a solution of 

S"^K-KS = VU. (B.4) 

Then 

V = -U(I + (HK)2)-1SV with 3 = S^e-*(^-^"') (B.5) 

solves the m x m discrete modified KdV equation (B.3). 
^^Note that V anticommutes with J, since U does, and we have = — V^. 
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Proof: From proposition 4.1 we recall that 

ZY = U(H-i-KHK)^iV, with H = S"e-'*(^+^"'-2i)J 

solves (B.l) if K is a solution of S~^K — KS = VU. Via (B.2), with the redefinitions 

S^-iSJ, V^iV, K^JK, 

(and keeping U) and using (4.2), this is translated into (B.4) and (B.5). □ 

The proof moreover gives a concrete map between a class of solutions of (B.l) and a class of 
solutions of (B.3). 

Appendix C: Equivalence of the fully discrete NLS equation with 
the corresponding equation of Ablowitz and Ladik in the scalar case 

After elimination of the terms cubic in q in favour of p, (6.7) reads 

i (<?+ -q) + {q^ - q)+ -{q- q~) + {pX -v)q+ + q (p+ - p~)* = o , 

p'^ — p = —eq'^q* . (CI) 
We can solve the second equation with 

p{x,t) = -e Y q{k,t)q{k-l,t)* x,teZ, (C.2) 

fc=— oo 

assuming the existence of this infinite series. Inserting this in the first equation, leads to 
i [q{x, t + l)- q{x, t)] +q{x+l,t+l) - q{x, t + 1) - [q{x, t) - q{x - 1, t)] 

X 

-eq{x,t) (q{k,t + iyq{k-l,t+l)-q{k-l,t)*q{k-2,t)^^ 

k=—oo 
x+1 

-e Y t + l)q{k-l,t + ly - q{k - 1, t) q{k - 2, i)*) q{x, t + 1) = . (C.3) 

fc=— oo 

Expressing this in terms of a new variable r related to q by g(x, t) = r(x, —t + l), and substituting 
t by —t in the resulting equation, in the scalar case we end up with equation (2.17) in [5] (with 
an obvious change in notation). The infinite sum in the last equation appears to be an awkward 
non-locality. But see [46] for how to achieve locality with the help of conservation laws. 

Appendix D: A non-local Hermitian reduction of the semi-dicrete 
NLS system 

Imposing the condition (see also [116]) 

q{x,t) = eq{-x,ty with e = ±1 (D.l) 

reduces the two equations (5.6) to the single equation 

i q{x, t) + q{x + l,t) - 2q{x, t) + q{x - l,t) - e q{x + l,t) q{-x, t^qix, t) 

-e q{x, t) q{-x, t^q^x - l,t) = . (D.2) 
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This is a non-local equation unless one imposes the restriction to functions q{x, t) that are even in 
X. Setting 

U = V^, U = eV^, S = {S-^)^, K^ = -K, = -K , (D.3) 
then H is given by with x replaced by —x, and proposition 5.1 implies that 

q = V^ (^{S^T e-^*'^(^^) - e K 5^ e^'^^^^^^K^ V , (D.4) 

with uj{S) defined in (5.10), solves (D.2) if K,K solve 

S-^K - K{S~^)^ = VV^ , S^K-KS = VVK (D.5) 

Appendix E: Lax pairs for the NLS systems 

(2.14) is the integrability condition of the linear equation 

d^r = (d(/.) ^' + (d^) A , (E.l) 

if A solves (cf. [1]) 

dA = (dA) A . (E.2) 

In the following we show how Lax pairs for the NLS systems originate from these equations, by 
using the respective bidifferential calculus. (E.2) is obviously satisfied if dA = dA = 0. These 
linear equations allow A to carry an arbitrary constant factor which turns out to play the role of 
a "spectral parameter" in the Lax pair. Although (E.l) is linear in this parameter, the resulting 
Lax pair exhibits a nonlinear parameter dependence. 

Continuous NLS system. (E.2) is trivially satisfied if dA = dA = 0. Using the bidifferential 
calculus that led to the continuous NLS system, this means A^; = A^ = [J, A] = 0. A particular 
solution is A = 2A / with a constant A. Then (E.l) becomes 

^^ = ]^J[J,^]^ + \[J,^, -i^t = J^x^ + 2\^^ (E.3) 
(where (j) = J(p)- Writing 

^ = ^e-^^'^-2i*^'^, (E.4) 
after some straight computations we obtain the linear system in the canonical form 

V'x = /:^, A = M^, (E.5) 

where 



£ = \J + U 



"A Im2 



M=U{U,-U^) + 2\\C = i( "^^ly ~ ,^^^1+ '^^^ 1 . (E.6) 
^ ' \ 2\q-q^ -2A2 1^^ + qq J ^ ^ 

This is a familar Lax pair for the matrix NLS system [3] . The integrability condition (zero curvature 
condition) 

Ct-M, + [C,M]=0 (E.7) 

indeed reproduces the NLS system (5.3). Whereas our linear system (E.l) is linear in the (spectral) 
parameter A, the above Lax pair has a nonlinear A-dependence. 
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Semi-discrete NLS system. (E.2) is trivially satisfied with A = 2aS ^, where a is a constant. 
Then (E.l) becomes 

= ^/^- + a[J,^], (E.8) 
-i4' = J (99+ - v?) ^' + 2 a - 

= J{U+ -U -U+U)"^ + 2a{'^+ - ^), (E.9) 



where we used the second of equations (3.19) in the last step. Multiplying from the right by 

2 



(l/2)(/ + J), and introducing = {I + J), we obtain 



'ip^=C'ilj, tp = M'il', (E.IO) 

where 

C = [{l + a)I -aJ]~^{I + l(), M = \J{U -U~ -UU~) + 2\a{C- I) . (E.ll) 

The integrability condition of this linear system, i.e. C = M'^ C — CM, is invariant under the 
transformations C = S+CS''^, M' = SMS''^ + 55"^ Setting 

(where x is the discrete variable), we recover the usual Lax pair of the AL system, 

\ q X w„,2 y ' V Xq - X {I - X ^) im^ + qq J 

Fully discrete NLS system. Choosing again A = 2a S~^, (E.l) yields 

=ZY^- + a[J,^r]^ -i (^'+ - *) = J(99+ - yj) ^+ + 2a (^ j - ^') . (E.14) 
Introducing again tj) = i^'~(/ + J) leads to 

V'+ = /:V, tp -ij:- = Mi: , (E.15) 

with 

C=\il+a)I -aJ]-^iI + U), M = ^—\J{^-^Z) + 2a{C-I)]. (E.16) 

i + 2a 

The integrability condition C — C- = L — C-M is invariant under the transformation C = 
S+CS-^, M' = I - S-{I - M)S"\ Choosing 



5 = A-(l-i-Ar'^ 0^ aLJ ' "=2^^ ^^) (^-1^) 
(where x,t £ Z), we obtain a Lax pair close to (E.13), 

c'=(^ {-^ ),M'=(^''" + i'"- + r p- ^ . - ^ - ) .(E.18) 

V Q X ^Im2 J V Ag„ - A (2 + 1 - A ^) 1^2 - P + P- y 
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Appendix F: Continuum limit 



Let us introduce lattice spacing parameters a,b in the basic relations that determine the bidiffer- 
ential calculus for the fully discrete NLS system, 

df = - [TS, f] Ci + ^[J, /] C2 , d/ = [T, /] Ci - - f] (2 , (F.l) 
a 2 a 

where now x E oL and t G 6Z. ^[TS, /] taken at a lattice point reads \[f{x + a,t + h) — 

f{x, t)] TS, which shows that an approach to the continuum can only work if we let first 6—7-0 and 
only afterwards a — >• 0. The first limit yields the bidifferential calculus for the semi-discrete NLS 
system and then the second that for the continuum NLS system (provided that / is differentiable) . 
In proposition 4.1 we only have to replace the expression for 3 by 

S = S^>(l-i^(S-i-I)jj (l-i^(I-S)J]-ij (F.2) 

fe->0 ga'/a g-it (S-i-I)J/a2 git (I-S)J/a2 _ gx/a ^-it (S+S"! -2I)J/a2 

In the last step we expressed S in terms of P via S = (I + aP)"^, which replaces (4.16). Renaming 
P to S, yields the expression for H in the continuous NLS case in proposition 4.1. 
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